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INTRODUCTION 

Since  the  invention  of  radar,  radar  operators  have  detected  and  tracked  targets  by  using  visual 
inputs  from  a  variety  of  displays.  Although  operators  can  perform  these  tasks  very  accurately,  they 
are  easily  saturated  and  quickly  become  fatigued.  Various  studies  have  shown  that  operators  can 
manually  track  only  two  to  four  targets.  To  correct  this  situation,  automatic  detection  and  tracking 
(ADI')  systems  were  attached  to  many  radars.  Undoubtedly,  as  digital  processing  increases  in  speed, 
and  hardware  decreases  in  cost  and  si/e,  ADT  systems  will  become  associated  with  almost  all  but  the 
simplest  radars. 

This  report  discusses  automatic  detection,  automatic  tracking,  and  sensor  integration  systems  for 
air-surveillance  radar.  Included  in  this  discussion  are  various  noncoherent  integrators  that  provide 
target  enhancement,  thresholding  techniques  for  reducing  false  alarms  and  target  suppression,  and 
algorithms  for  estimating  target  position  and  resolving  targets.  Then,  an  overview  is  given  of  the 
entire  tracking  system  followed  by  a  discussion  of  its  various  components  such  as  tracking  filter, 
maneuver-following  logic,  track  initiation,  and  correlation  logic.  Next,  multiscan  approaches  to 
automatic  tracking  such  as  maximum  likelihood  are  discussed.  Finally,  the  report  concludes  with  a 
discussion  of  sensor  integration  and  radar  netting,  including  colocated  and  multisite  systems. 

AUTOMATIC  DETECTION 

In  the  ld40s,  Marcum  [1]  introduced  to  the  radar  community  the  statistical  framework  necessary 
for  the  development  of  automatic  detection.  Later,  Swerling  [2]  extended  the  work  to  fluctuating 
targets.  They  investigated  many  of  the  statistical  problems  associated  with  the  noncoherent  detection 
of  targets  in  Rayleigh  noise.  The  most  important  result  of  the  investigation  was  the  generation  of 
curves  of  probability  of  detection  (Pn)  vs  signal-to-noise  ratio  (S/N)  for  a  detector  that  sums  N 
enveloped-detected  samples  (cither  linear  or  square  law)  under  the  assumption  of  equal  signal  ampli¬ 
tudes.  However,  in  a  search  radar,  as  the  beam  sweeps  over  the  target,  the  returned  signal  amplitude 
is  modulated  by  the  antenna  pattern.  Many  authors  investigated  various  detectors  (weightings),  com¬ 
paring  detection  performance  and  angular  estimation  results  to  the  optimal  values. 

In  the  original  work  on  these  detectors,  the  environment  was  assumed  to  be  known  and  homo¬ 
geneous  so  that  fixed  thresholds  could  be  used.  However,  a  realistic  environment  (e.g.,  containing 
land,  sea,  and  rain)  will  cause  an  exorbitant  number  of  false  alarms  for  a  fixed  threshold  system  that 
does  not  use  excellent  coherent  processing;  therefore,  three  main  approaches,  adaptive  thresholding, 
nonparametric  detectors,  and  clutter  maps  were  used  to  solve  the  false-alarm  problem.  Adaptive 
thresholding  and  nonparametric  detectors  arc  based  on  the  assumption  that  homogeneity  exists  in  a 
small  region  about  th?  range  cell  that  is  being  tested.  The  adaptive  thresholding  method  assumes  that 
the  noise  density  is  known  except  for  a  few  unknown  pnmmeters  (c.g.,  the  mean  and  variance).  The 
surrounding  reference  cells  are  then  used  to  estimate  the  unknown  parameters,  and  a  threshold  is 
obtained  based  on  the  estimated  density.  Nonparametric  detectors  obtain  a  constant  false-alarm  rate 
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(t  FAR)  by  ranking  (i.e.,  ordering  the  samples  trom  smallest  to  largest)  the  test  samples,  usually  with 
the  lelerence  cells.  Under  the  hypothesis  that  al  the  samples  (test  anil  reference)  are  independent 
samples  trom  an  unknown  density  function,  the  test  sample  has  a  uniform  density  function;  conse¬ 
quently  a  threshold  that  yields  CTAR  can  he  set.  Clatter  maps  store  an  average  background  level  for 
each  range-azimuth  cell.  A  target  is  then  declared  in  a  n:nge-azimuth  cell  if  the  new  value  exceeds 
the  average  background  level  by  a  specified  amount. 


Optima!  Detector 

The  radar-detection  problem  is  a  binary  hypothesis  testing  problem 


No  target  present 


H | :  Target  present. 

Several  criteria  (i.e.,  definitions  of  optimality)  can  be  used  to  solve  this  problem,  but  the  most 
appropriate  for  radar  is  the  Neyman-Pearson  (3).  This  criterion  maximizes  the  Pt)  for  a  given  proba¬ 
bility  of  false  alarm  P!u  by  comparing  the  likelihood  of  ratio  L  to  an  appropriate  threshold  T.  A  tar¬ 
get  is  declared  present  if 


Uv,. 


P<*\ . Xn  \  H\) 

pU\ . \H{)) 


(1) 


where  p (.v , . v„  |  //,)  and  (.v ! . x„  j  Hf) )  are  the  joint  densities  of  the  n  samples  .t,  under  the 

conditions  of  target  presence  and  target  absence  respectively.  For  a  linear  envelope  detector  the  sam¬ 
ples  have  a  Rayleigh  density  under  H{)  and  a  Ricean  density  under  H  |,  and  the  likelihood  detector 
reduces  to 


n 


L  A, 


AjXj 

cr 


>  T. 


(2) 


where  /()  is  the  Bessel  function  of  zero  order,  <r  is  the  noise  power,  and  At  is  the  target  amplitude  on 
the  /th  pulse  and  is  proportional  to  the  power  antenna  pattern.  For  small  signals  (A,  «  a),  the 
detector  reduces  to  the  square-law  detector; 


E  Arx  'r  >  t. 


(3) 


For  Aj  -  A,  this  detector  and  the  linear  detector  were  studied  first  by  Marcum  [l|  and  in  the 
succeeding  years  by  numerous  other  people.  The  most  important  facts  concerning  these  detectors  are 
as  follows: 


•  The  detection  performances  of  the  linear  and  square-law  detectors  are  similar.  The  difference 
is  less  than  0.2  dB  over  wide  ranges  of  P  n •  P fa '  and  n  . 


:>>  N 


Since  the  siunai  return  of  a  scanning  '  ;■  ;r  tulated  by  the  ,n!,  nna  put •>' 

only  tt.K4  of  the  pulses  between  the  halt-power  points  should  be  integrated,  and  th 


iua.\.:v,'/.e 
e  antenna 
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beam  shape  factor  (ANSI*')  is  l.f>  dll  |4|.  The  A H S I •  is  the  number  by  which  the  midbeam  S/N  must 
be  reduced  so  that  the  detection  curves  generated  Cor  equal  signal  amplitudes  can  be  used  for  the 
scanning  radar. 

•  Figure  1  shows  that  the  collapsing  loss  lor  the  linear  detector  can  be  several  dB  greater  than 
the  loss  for  a  square-law  detector  |5|.  The  collapsing  loss  is  the  additional  signal  required  to  maintain 
the  same  and  /’M  when  unwanted  noise  samples  along  with  the  desired  signal-plus-noisc  samples 
are  integrated.  The  number  of  signal  samples  integrated  is  N ,  the  number  of  extraneous  noise  sam¬ 
ples  integrated  is  A/,  and  the  collapsing  ratio  p  -=  ( N  +  M)/N . 


Fig.  I  —  Collapsing  loss  vs  collapsing  ratio  for  a  probability  of  false  alarm  of  10  6 
and  a  probability  of  detection  of  0.5,  (copyright  1972,  IEEE).  From  G.V.  Trunk 
15] 


Most  automatic  detectors  are  required  not  only  to  detect  targets  but  to  make  angular  estimates  of 
the  azimuth  position  of  the  target.  Swerling  [6]  calculated  the  standard  deviation  of  the  optimal  esti¬ 
mate  by  using  the  Cramer-Rao  lower  bound.  The  results  are  shown  in  Fig.  2,  where  a  normalized 
standard  deviation  is  plotted  against  the  midbeam  S/N.  This  result  holds  for  a  moderate  or  large 
number  of  pulses  integrated,  and  the  optimal  estimate  involves  finding  the  location  where  the  correla¬ 
tion  of  the  returned  signal  and  the  derivative  of  the  antenna  pattern  is  zero.  Although  this  estimate  is 
rarely  implemented,  its  performance  is  approached  by  simple  estimates. 


_2  0  2  4  6  8  10  12 

S/N  (dB) 

Fig.  2  -  Comparison  of  angular  estimates  with  the  Cramer-Rao  lower  bound. 
a  is  the  standard  deviation  of  the  estimation  error,  and  N  is  the  number  of 
pulses  within  the  3  dB  bcamwidth,  which  is  2  (3/0.85.  The  S/N  is  the  value 
at  the  center  of  the  beam,  (copyright  1970,  IEEE).  From  G.V,  Trunk  |llj 
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Practical  Detectors 

Many  different  detectors  are  used  to  accumulate  the  radar  returns  as  a  radar  sweeps  by  a  target. 
Figure  3  shows  a  few  of  the  most  common  detectors  [7].  Though  they  arc  shown  as  being  imple¬ 
mented  with  shift  registers,  normally  they  arc  implemented  with  random  access  memory.  The  input 
to  these  detectors  can  be  either  linear,  square-law,  or  log  video.  Since  linear  is  probably  the  most 
commonly  used,  the  advantages  and  disadvantages  of  the  various  detectors  will  be  stated  for  this 
video. 


INTEGRATORS 


(a)  Moving  window 


(b)  Feedback  integrator 
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(c)  Two-Pole  filter 


(d)  Binary  integrator 


INTEGDATE  I 
K  PULSES,  - 
THEN  DUMP  I 


BINARY  | 
INTEGRATOR  ‘ 


(e)  Batch  processor 

Fig.  3  —  Block  diagrams  of  various  detec¬ 
tors.  The  ielter  C  indicates  a  comparison 
and  t  indicates  a  delay.  From  G.V.  Trunk 
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Moving  Window 

The  moving  window  in  Fig.  3(a)  performs  a  running  sum  of  n  pulses  in  each  range  cell; 

S,  =  Si  _  I  +  X;  -  Xi_„,  (4) 

or  equivalently. 

Si  =2:  (5) 

y=0 

where  5,  is  the  sum  on  the  / th  pulse.  The  performance  of  this  detector  is  only  0.5  dB  worse  than  the 
optimal  detector  given  by  Eq.  (3)  [8],  The  detection  performance  can  be  obtained  by  using  an  ABSF 
of  1.6  dB  and  the  detection  curves  in  Ref.  1.  The  angular  estimate  that  is  obtained  by  either  taking 

the  maximum  value  of  the  running  sum  or  by  taking  the  midpoint  between  the  first  and  last  crossing 

of  the  detection  threshold  has  a  bias  of  n/2  pulses  that  can  be  easily  corrected.  The  standard  devia¬ 
tion  of  the  estimation  error  of  both  estimators  is  about  20%  higher  than  that  for  the  optimal  estimate. 

The  major  disadvantage  of  this  detector  is  its  susceptibility  to  interference;  that  is,  one  large  sample 

from  interference  can  cause  a  detection.  A  minor  disadvantage  is  that  the  last  n  pulses  for  each  range 
cell  must  be  saved,  and  therefore  a  large  storage  is  required  when  a  large  number  of  pulses  are 
integrated.  However,  because  of  the  availability  of  large  memories  of  reduced  size  and  cost,  this  is  a 
minor  problem. 

The  detection  performance  discussed  is  based  on  the  assumption  that  the  target  is  centered  in  the 
moving  window.  In  a  real  situation  the  radar  scans  over  the  target  and  decisions  that  are  highly 
correlated  are  made  at  every  pulse.  Hansen  [9]  analyzed  this  situation  and  calculated  the  detection 
thresholds  shown  in  Fig.  4,  the  detection  performance  shown  in  Fig.  5,  and  the  angular  accuracy 
shown  in  Fig.  6.  Comparing  Hansen’s  calculation  with  the  single  point  calculation,  one  concludes 
that  1  dB  of  improvement  is  obtained  by  making  a  decision  at  every  pulse.  Figure  6  shows  that  the 
angular  accuracy  of  the  beam-splitting  procedure  is  about  20%  greater  th„n  the  optimal  estimate.  For 
large  signal-to-noise  ratios,  the  accuracy  (rms  error)  of  the  beam-splitting  and  maximum  return  pro¬ 
cedures  will  be  limited  by  the  pulse  spacing  [10]  and  will  approach 

a  (0)  =  AO/'Tm  ,  (6) 

where  AO  is  the  angular  rotation  between  transmitted  pulses.  Consequently,  if  the  number  of  pulses 
per  beamwidth  is  small,  the  angular  accuracy  will  be  poor.  For  instance,  if  pulses  are  separated  by 
0.5  beamwidth,  a(0)  is  bounded  by  0.14  beamwidth.  However,  improved  accuracy  can  be  obtained 
by  using  the  amplitudes  of  the  radar  returns.  If  one  assumes  a  Gaussian-shaped  antenna  pattern,  an 
accurate  estimate  of  the  target  angle  is  given  by 

0  =  0 ,  +  InMj/^t)  (7) 

2  2 a  AO 

where 

a  =  1 ,386/[beamwidthJ2 ,  (8) 

and  A\  and  A2  are  the  two  largest  amplitudes  of  the  returned  samples  and  occur  at  angles  6\  and 
02  =  0)  +  Ad  respectively.  Since  the  estimate  should  lie  between  0(  and  02  but  Eq.  (7)  will  not 
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Fig.  4  —  Single-sweep  false-alarm  probability  Pfa  vs  threshold  for 
moving  window,  (copyright  1970,  IEEE).  From  V.G.  Hansen  (9) 


SINGLE-PULSE  SIGNAL-TO-NOISE  RATIO  (dB) 


Fig.  5  —  Detection  performance  of  the  ana¬ 
log  moving  window  detector  for  the  no  fad¬ 
ing  case,  (copyright  1970,  IEEE).  From 
V.G.  Hansen  19] 


Fig.  b  —  Angular  accuracy  obtained  with 
beam-splitting  estimation  procedure  for  the  no- 
ladtng  ease.  Broken  line  curves  are  lower 
hounds  derived  by  Swerling  (61,  and  points 
shown  are  simulation  results,  (copyright  1970. 
IEEE).  From  V.G.  Hansen  |9| 
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always  yield  such  as  estimate,  0  should  be  set  equal  to  0,  if  0  <  0,  and  0  should  be  set  equal  to  6 2  if 
6  >  0;.  Figure  7  shows  the  accuracy  of  this  estimator  for  the  case  of  n  =  2  pulses  per  beamwidth. 
This  estimation  procedure  can  also  be  used  to  estimate  the  elevation  angle  of  a  target  in  multibeam 
systems. 


S'N  idB) 


Fig  7  —  Angular  accuracy  using  two-pulse  estimates 


Feedback  Integrator 

The  amount  of  storage  required  can  be  reduced  significantly  by  using  a  feedback  integrator  as 
shown  in  Fig.  3(b) 


j,  =  K  Sk~)  +  x,  ■  (9) 

For  a  feedback  vaiue  of  K .  the  effective  number  of  pulse  integrated  M  is  M  —  1/(1  -  K)  and  for 
optimal  (maximum  Pr) )  performance  M  =  0.63  N,  where  .V  is  the  number  of  pulses  between  the 
3  dB  antenna  beamwidth  |llj  The  detection  performance  is  given  by  the  detection  curves  for  M 
pulses  with  ABSF  =  1.6  dB  1 1  I).  Although  the  feedback  integrator  applies  an  exponential  weighting 
into  the  past,  its  detection  performance  in  only  I  dB  less  than  the  optimal  integrator  |8).  However, 
difficulties  are  encountered  when  using  the  feedback  integrator  to  estimate  the  azimuth  position  111] 
The  threshold  crossing  procedure  yields  estimates  only  209c  greater  than  the  lower  bound,  but  the 
bias  is  a  function  of  the  S.  N  and  must  be  estimated  On  the  other  hand,  the  maximum  value,  though 
it  has  a  constant  bias,  has  estimates  that  are  100%  greater  than  those  of  the  lower  bound.  Further¬ 
more,  the  exponential  weighting  function  essentially  destroys  the  radar  sidelobes.  Because  of  these 
problems,  the  feedback  integrator  is  limited  in  usefulness. 
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I  lie  two-polo  t titer  in  Fig.  3(c)  requires  the  storage  of  an  intermediate  calculation  besides  the 
integrated  output  and  is  described  mathematically  by 

v,  =  x,  —  k2  z,-_i  (10) 

and 

—  y ,  _  i  +  &  1  Z/  _  i ,  (11) 

where  xt  is  the  input,  y-  is  the  intermediate  calculation,  z,  is  the  output,  and  kt  and  k2  are  the  two 
feedback  values.  The  values  [12,13]  that  maximize  Pn  are  given  by 

*i=2  exp  ( —  1  —  T2)  cos  ( <j>dT )  (12) 

and 

k2  =  exp  (-2  -  f2) ,  (13) 

where  f  =  0.63,  Nojjt  —  2.2,  and  N  is  the  number  of  pulses  between  the  3  dB  points  of  the 
antenna.  With  this  rather  simple  device,  a  weighting  pattern  similar  to  the  antenna  pattern  can  be 
obtained.  The  detection  performance  is  within  0.15  dB  of  the  optimal  detector,  and  its  angular  esti¬ 
mates,  as  shown  in  Fig.  2,  are  about  20%  greater  than  the  Cramer-Rao  lower  bound.  If  the  desired 
number  of  pulses  integrated  is  changed  because  the  rotation  of  the  radar  is  changed  or  because 
another  radar  is  used,  it  is  necessary  only  to  change  the  feedback  values  and  k2.  The  problems  of 
this  detector  arc  that  (a)  it  has  rather  high  detector  sidclobes,  15  to  20  dB  and  (b)  it  is  extremely  sen¬ 
sitive  to  interference  (i.e.,  the  filter  has  a  high  gain  that  results  in  a  large  output  for  a  single  sample 
that  has  a  high  value). 

Binary  Integrator 

The  binary  integrator  is  also  known  as  the  dual-threshold  detector,  Af-out-of-A/  detector,  or  rank 
detector  (see  section  on  Nonparametric  Detectors).  Numerous  individuals  have  studied  it  [14-18],  As 
shown  in  Fig.  3(d),  the  input  samples  are  quantized  to  0  or  1  depending  on  whether  or  not  they  are 
less  than  a  threshold  7",.  The  las'  N  zeroes  and  ones  are  summed  and  compared  to  a  second  thres¬ 
hold  7\  =  M .  For  large  N .  the  detection  performance  of  this  detector  is  approximately  2  dB  less 
than  that  of  the  moving-window  integrator  because  of  the  hard  limiting  of  the  data,  and  the  angular 
estimation  error  is  about  25%  greater  than  the  Cramer-Rao  lower  bound.  Schwartz  [16]  showed  that 
within  0.2  dB,  the  optimal  value  of  M  for  maximum  PD  is  given  by 

M  =  1.5  (14) 


when  10'" 10  <  Pfu  <  10“5  and  0.5  <  PD  <  0.9.  Dillard  [18],  as  shown  in  Fig.  8,  calculated  the 
optimal  value  of  PN ,  and  the  probability  of  exceeding  Tx  when  only  noise  is  present.  Then,  the 
threshold  7’,  is  found  from 


OO  v 

*  I,  7 ' 


-,r:/2  a2 


dx  , 


(15) 
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Fig.  8  —  Optimum  values  of  /\  as  a  function  of  the  sample  size,  n ;  probability  of 
false  alarm,  or;  Rice  distribution,  0=1.  (copyright  1967,  IEEE).  From  G.M.  Dil¬ 
lard  1 18) 


which  yields 

T\  =  o  ( —  2  In  .  (16) 

Figures  9  and  10  give  a  comparison  of  the  optimal  (best  value  of  M )  binary  integrator  with  various 
other  procedures. 

The  binary  integrator  is  used  in  many  radars  because  (a)  it  is  easily  implemented,  (b)  it  ignores 
interference  spikes  that  cause  trouble  with  integrators  directly  using  signal  amplitude,  and  (c)  it  works 
extremely  well  when  the  noise  has  a  non-Rayleigh  density  [19],  Figure  1 1  shows  a  comparison  of  the 
optimal  binary  integrator  (three-out-of-three,  M-out-of-A^),  another  binary  integration  (two-out-of- 
threc).  and  the  moving  window  detector  in  log-normal  interference  (an  example  of  a  non-Rayleigh 
density).  The  optimal  binary  integrator  is  much  better  than  the  moving  window  integrator.  The 
optimal  values  for  log-normal  interference  were  calculated  by  Schleher  [19]  and  are  m  =  3,  8,  and 
25.  lor  n  =  3,  10,  and  30  respectively. 

A  modified  version  of  binary  integration  is  used  sometimes  when  there  are  a  large  number  of 
pulses.  It  also  h',s  flexibility  to  integrate  a  different  number  of  pulses.  The  modified  binary  moving 
window  (MBMW)  differs  from  the  ordinary  binary  moving  window  (OBMW)  because  of  the  intro¬ 
duction  of  a  third  threshold.  When  the  second  threshold  is  reached,  one  counts  the  number  of  con¬ 
secutive  pulses  for  which  the  second  threshold  is  exceeded.  When  this  number  equals  the  third  thres¬ 
hold,  a  target  is  declared.  The  performance  of  the  MBMW  and  a  comparison  to  the  OBMW  are 
shown  in  Table  1  [20],  The  larger  the  value  of  N ,  the  larger  is  the  difference  in  dB  performance 
between  the  MBMW  and  OBMW  detectors. 
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S/N  (dB) 


Fig.  9  —  Comparison  of  binary  iniegrators  (M-out-of-iV)  with 
other  integration  methods  ( P,„  =  10  “l0.  P»  =  0.5),  (copyright 
1956.  IEEE).  From  M.  Schwartz  (161 


Fig.  10  —  Comparison  of  binary  integrator  (M-out-of-(V)  with 
other  integration  methods  (P,„  =  10  Pn  =  .90).  From 
M.  Schwartz  1 16J 
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Fig.  II  -  Comparison  of  various  detectors  in  log-normal  interfer¬ 
ence  (,V  =  3,  Pfu  -  10  h),  (copyright  1975,  lEEF.).  After  D.C. 
Schleher  |19| 


Table  1  —  Comparisons  of  the  Detection  Performance  Between 
MBMW  and  OBMW  at  A  =  8,  16,  24, 

32  with  P,  =  10~\  Pn  =  0.9.  After  Mao  [20] 


Number  of 

Type  of 

Detection 

T, 

S/N(dB) 

Pulses 

Detector 

Criteria 

r  N 

1  1 

ll 

00 

MBMW 

4/7  +  1 

0.018 

3.22 

8.28 

OBMW 

5/8 

0.032 

3.015 

8.13 

A  =  16 

MBMW 

4/7  +  5 

0.054 

3.81 

7.54 

MBMW 

4.  3/7  +  7 

0.060 

2.77 

7.01 

OBMW 

8/16 

0.064 

2.74 

6.21 

MBMW 

4/7  +  6 

0.071 

2.70 

7.34 

A  =  24 

MBMW 

4,  3/7  +  7 

0.060 

2.77 

6.98 

OBMW 

10/24 

0.072 

2.69 

5.18 

MBMW 

4/7  +  7 

0.091 

2.59 

7.32 

A  =  32 

MBMW 

4,  3/7  +  7 

0.060 

2.77 

6.98 

OBMW 

13/32 

0.095 

2.57 

4.53 

Batch  Processor 

The  batch  processor  is  very  useful  when  a  large  number  of  pulses  are  in  the  3  dB  beamwidth. 
If  KN  pulses  arc  in  the  3  dB  beamwidth  as  shown  in  Fig.  3,  K  pulses  are  summed  (batched)  and 
either  a  0  or  1  is  declared  depending  on  whether  or  not  the  batch  is  less  than  a  threshold  7,.  The  last 
A  zeros  and  ones  are  summed  and  compared  to  a  second  threshold  M .  An  alternative  version  of  this 
detector  is  to  put  the  batches  through  a  moving  window  detector. 
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The  batch  processor,  like  the  binary  integrator,  is  easily  implemented,  it  ignores  interference 
spikes,  and  works  extremely  well  when  the  noise  has  a  non-Rayleigh  density.  Furthermore,  when 
compared  with  the  binary  integrator,  the  batch  processor  requires  less  storage,  detects  better,  and  esti¬ 
mates  angles  more  accurately.  For  instance,  if  there  were  80  pulses  on  target,  one  could  batch  16 
pulses,  quantize  this  result  to  a  0  or  1,  and  declare  a  target  with  a  3-out-of-5  (or  2-out-of-5)  binary 
integrator.  With  an  8  bit  analog-to-digital  converter,  the  storage  requirement  per  range  cell  is  17  bits 
(12  bits  for  batch  and  5  for  binary  integrator)  for  the  batch  processor  as  opposed  to  80  bits  for  the 
binary  integrator  and  640  bits  for  the  moving  window.  The  detection  performance  of  the  batch  pro¬ 
cessor  tor  a  large  number  of  pulses  integrated  is  approximately  0.5  dB  worse  than  the  moving  win¬ 
dow. 


The  batch  processor  has  been  successfully  implemented  by  the  Applied  Physics  Laboratory 
(APL)  (21)  of  Johns  Hopkins  University.  To  obtain  an  accurate  azimuth  estimate  approximately  20 7c 
greater  than  the  lower  bound,  APL  uses 


LB,  ' 


(17) 


where  Bt  is  the  batch  amplitude  and  0,  is  the  azimuth  angle  that  corresponds  to  the  center  of  the 
batch . 


False  Alarm  Control 

In  the  presence  of  clutter,  if  fixed  thresholds  are  used  with  the  previously  discussed  integrators, 
an  enormous  number  of  detections  will  occur  and  will  saturate  and  disrupt  the  Hacking  computer  asso¬ 
ciated  with  the  radar  system.  Four  important  facts  should  be  noted: 

•  It  makes  little  sense  to  have  an  automatic  detection  system  without  an  associated  tracking  sys¬ 
tem; 

•  the  of  the  detector  should  be  as  high  as  possible  without  saturating  the  tracking  computer; 

•  Random  false  alarms  and  unwanted  targets  (e  g.,  stationary  targets)  arc  not  a  problem  if  they 
are  removed  bv  the  tracking  computer;  and 

•  scan-to-scan  processing  is  the  only  way  to  remove  stationary  point  clutter  or  moving  target 
indicator  (MT1)  clutter  residues. 

One  can  limit  the  number  of  false  alarms  with  a  fixed  threshold  system  by  setting  a  very  high  thres¬ 
hold.  However,  this  would  reduce  target  sensitivity  in  regions  ot  low  noise  (clutter)  return.  To 
reduce  the  false-alarrn  problem,  three  main  approaches  have  been  used:  adaptive  threshold,  non- 
parametric  detectors,  and  clutter  maps.  Adaptive  thresholding  and  nonparametric  detectors  assume 
that  the  samples  in  the  range  cells  surrounding  the  test  cell  (called  reference  cells)  arc  independent 
and  identically  distributed  Usually  it  is  assumed  that  the  time  samples  are  independent.  Both  kinds 
ol  deicctois  test  whelncr  the  test  cell  has  a  return  sufficiently  larger  than  the  reference  cells.  Clutter 
maps  allow  variation  in  space  bui  require  stationarity  over  several  scans.  Clutter  maps  store  an  aver¬ 
age  background  level  for  each  range-u/imuth  cell.  A  target  is  then  declared  in  a  range-azimuth  cell  if 
the  new  value  exceeds  the  average  background  level  by  a  specified  amount. 
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1  be  basic  assumption  oi  the  adaptive  thresholding  technique  is  that  the  probability  density  of  the 
noise  is  known  except  lor  a  lew  unknown  parameters.  The  surrounding  reference  cells  are  used  to 
estimate  the  unknown  parameters,  anti  a  threshold  based  on  the  estimated  parameters  is  obtained.  The 
simplest  adaptive  detector,  shown  in  Fig.  12,  is  the  cell-averaging  CFAR  investigated  by  Finn  and 
Johnson  122],  II  the  noise  has  a  Rayleigh  density,  />(.v)  =  x  exp  ( —  x~ /2cr)/cr ,  only  the  parameter  a 
needs  to  be  estimated  and  the  threshold  is  ol  the  form  T  —  K^Zxj  =  Kn'^ir/2a.  Flowever,  since  7  is 
set  by  an  estimate  a,  it  has  some  error  anti  must  be  slightly  larger  than  the  threshold  one  would  use  if 
|T  were  known  exactly  a  priori.  The  raised  threshold  causes  a  loss  in  target  sensitivity  and  is  referred 
to  as  a  CFAR  loss.  This  loss  has  been  calculated  [23|  and  is  summarized  in  Table  2.  Table  2  shows 
that  tor  a  small  number  of  reference  cells  the  loss  is  large  because  of  the  poor  estimate  of  a. 
Consequently,  a  large  number  of  reference  cells  would  be  preferred.  However,  if  this  occurs,  the 
homogeneity  assumption  (i.e. .  all  the  reference  cells  are  statistically  similar)  might  be  violated.  A 
good  rule  of  thumb  is  to  use  enough  reference  cells  so  that  the  CFAR  loss  is  below  1  dB  and  at  the 
same  time  not  let  the  reference  cells  extend  beyond  1  mi  on  either  side  of  the  test  cell.  For  a  particu¬ 
lar  radar  this  might  not  be  feasible. 


II  there  is  uncertainty  about  whether  or  not  the  noise  is  Rayleigh  distributed,  it  is  better  to  thres¬ 
hold  individual  pulses  and  use  a  binary  integrator  as  shown  in  Fig.  13.  This  detector  is  tolerant  to 
variations  in  the  noise  density  because  by  setting  K  to  yield  a  “1”  with  probability  0.1,  a 
PUt  ~  10  fl  can  be  obtained  by  using  a  7-out-of-9  detector.  Although  noise  may  be  non-Rayleigh,  it 
will  probably  be  very  Rayleigh-like  out  to  the  tenth  percentile.  Furthermore,  one  can  use  feedback 
based  on  several  scans  of  data  to  control  K  in  order  to  maintain  a  desired  either  on  a  scan  or  sec¬ 
tor  basis,  thus  demonstrating  the  general  rule  that  to  maintain  a  low  Pja  in  various  environments, 
adaptive  thresholding  should  be  placed  in  front  of  the  integrator. 


It  the  noise  power  varies  from  pulse-to-pulse  (as  it  would  in  jamming  when  frequency  agility  is 
employed),  one  must  CFAR  each  pulse  and  then  integrate.  While  the  binary  integrator  performs  this 
type  of  CFAR  action,  analysis  [24,25]  has  shown  that  the  ratio  detector  shown  in  Fig.  14  is  a  better 
detector.  The  ratio  detector  sums  S/N  ratios  and  is  specified  by 


n 


£ 

/  =  i 


xi2  O') 


I  m 

~  E  IvU  +  1  +  k)  +  X,2(J  -  1  -  *)] 
2m  k  =  t 


(18) 


where  x,{j)  is  the  / th  envelope-detected  pulse  in  the  y'th  range  cell  and  the  number  of  reference  cells 
is  2  m.  The  denominator  is  the  maximum  likelihood  estimate  of  <j,2,  the  noise  power  per  pulse.  It 
will  detect  targets  even  though  only  a  few  returned  pulses  have  a  high  S/N  ratio.  Unfortunately,  this 
will  also  cause  the  ratio  detector  to  declare  false  alarms  in  the  presence  of  narrow-pulse  interference. 
To  reduce  the  number  of  false  alarms  when  narrow-pulse  interference  is  present,  the  individual  power 
ratios  can  be  soft-limited  [25|  to  a  small  enough  value  so  that  interference  will  only  cause  a  few  false 
alarms.  Figures  15  and  16  show  a  comparison  of  the  ratio  detector  with  other  commonly  used  detec¬ 
tors.  Figure  17  shows  a  typical  performance  in  sidclobc  jamming  when  the  jamming  level  varies  by 
20  dB  per  pulse.  By  employing  a  second  test  to  identify  the  presence  of  narrow-pulse  interference,  a 
detection  performance  approximately  halfway  between  the  limiting  and  nonlimiting  ratio  detectors  can 
be  obtained. 


If  the  noise  samples  are  dependent  in  time  or  have  a  non-Rayleigh  density,  such  as  the  chi- 
square  density  or  log-normal  density,  it  is  necessary  to  estimate  two  parameters  and  the  adaptive 
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Fig.  12  —  Cell-averaging.  From  CFAR  (7) 


Table  2  —  CFAR  Loss  for  Pfu  =  10-6  and  PD  =  0.9. 
After  Mitchell  and  Walker  [23] 


Number  of 
Pulses 
Integrated 

Loss  for  Various  Numbers 
of  Reference  Cells  (dB) 

1 

2 

3 

5 

10 

00 

1 

— 

— 

15.3 

7.7 

3.5 

0 

3 

7.8 

5.1 

3.1 

1.4 

0 

10 

6.3 

3.3 

2.2 

1.3 

0.7 

0 

30 

3.6 

2.0 

1.4 

1.0 

0.5 

0 

100 

2.4 

1.4 

1.0 

0.3 

0 
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Fig  13  —  Implementation  of  binary  integrator.  From  G.V.  Trunk  (7) 
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Fig  14  —  Ratio  detector 
From  G.V.  Trunk  (71 
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Fig  15  —  Curves  of  probability 
of  detection  vs  signal-to-noise 
ratio  per  pulse  for  the  cell¬ 
averaging  CFAR.  rttio  detectors, 
log  integrator,  and  binary 
integrator:  nonfluctuating  target 
and  probability  of  false  alarm 
=  lO"6  From  G.V.  Trunk  [25 ) 
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Fig  16  —  Curves  of  probabilit>  of  detection  vs  signal-lo-noise  ratio  for  the  cell¬ 
averaging  CFAR,  ratio  detectors,  log  integrator,  and  binary  integrator:  Rayleigh 
fluctuating  target  and  probability  of  false  alarm  -  10  "6.  From  G.V  Trunk  (25] 


MAXIMUM  JAMMING /NOISE  ■  20  dB 
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Fig.  17  —  Curves  of  probability  of  detection  vs  signal-to-noisc  ratio  for  the  celi¬ 
ac  craging  CFAR.  ratio  detectors,  log  integrator,  and  binary  integrator  Rayleigh 
fluctuations,  probability  of  false  alarm  =  10  *  and  maximum  jamming-io-noise  ratio 
=  20  dB.  From  G.V.  Trunk  (25] 


detector  is  more  complicated.  Usually  several  pulses  are  integrated  so  that  one  can  assume  that  the 
integrated  output  has  a  Gaussian  probability  density  .  Then  the  two  parameters  that  must  be  estimated 
are  the  mean  and  variance,  and  a  threshold  of  the  form  T  =  g  +  A  a  is  used.  Though  the  mean  is 
easilv  obtained  in  hardware,  the  usual  estimate  of  the  standard  deviation. 


=  77  -  <)2 


where 


S  2  * 


is  more  difficult  to  implement.  Consequently,  the  mean  deviate  defined  by 


NRL  REPORT  9110 


<T  =  A  S  1  .V; 


(21) 


is  sometimes  used  because  it  is  easy  to  implement.  Nothing  can  be  done  to  the  binary  integrator  to 
yield  a  low  Pfll  it  the  noise  samples  are  eorrelated.  Thus,  it  should  not  be  used  in  this  situation. 
However,  it  the  correlation  time  is  less  than  a  batching  interval,  the  batch  processor  will  yield  a  low 
P,Ll  without  modifications. 

Target  Suppression 

Target  suppression  is  the  loss  in  detectability  that  is  caused  by  other  targets  or  clutter  residues  in 
the  reference  cells.  Basically  there  are  two  approaches  for  solving  this  problem:  (1)  remove  large 
returns  Irom  the  calculation  of  the  threshold  [26,27);  (2)  diminish  the  effects  of  large  returns  by  either 
limiting  or  using  log  video.  The  technique  that  should  be  used  is  a  function  of  the  particular  radar 
system  and  its  environment. 

Rickard  and  Dillard  [27 ]  proposed  a  class  of  detectors  DK  where  the  K  largest  samples  are  cen¬ 
sored  from  the  reference  cells.  A  comparison  of  D0  (no  censoring)  with  D  |  and  D2  for  a  Swerling  2 
target  is  shown  in  Fig.  18.  where  N  is  the  number  of  reference  cells,  (3  is  the  ratio  of  the  power  of 
the  interfering  target  to  the  target  in  the  test  cell,  and  the  bracketed  pair  (m  ,n)  indicates  the  Swerling 
models  of  the  target  and  the  interfering  target  respectively.  Figure  18  shows  that  when  one  has  an 
interfering  target  the  Pp  does  not  approach  1  as  S/N  increases.  Another  approach  [26],  which  cen¬ 
sors  samples  in  the  reference  cell  if  they  exceed  a  threshold,  is  briefly  discussed  in  the  section  about 
nonparametric  detectors. 
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Fig.  18  —  Detection  probability  vs  SNR  for  Swerling  case  2  primary  tar¬ 
get  (copyright  1977,  IEEE).  From  J.T.  Rickard  and  G.M.  Dillard  [27| 


The  use  of  log  video  is  an  alternate  approach  for  interfering  targets.  Basically  by  taking  the  log, 
large  samples  in  the  reference  cells  will  have  little  effect  on  the  threshold.  The  loss  associated  with 
the  use  of  log  video  is  0.5  dB  for  10  pulses  integrated  and  1.0  dB  for  100  pulses  integrated  [28J. 
Figure  19  shows  an  implementation  of  the  log-CFAR  [29],  To  maintain  the  same  CFAR  loss,  the 
number  of  reference  cell  A/log  for  the  log-CFAR  should  equal 

A/|„g  =  1.65  MUn  —  0.65,  (22) 
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Fig.  19  —  Block  diagram  of  cell  averaging  LOG/CFAR  receiver,  (copy¬ 
right  1972,  IEEE).  From  V.G.  Hansen  and  J.R.  Ward  |29j 

where  MUn  is  the  number  of  reference  cells  for  linear  video.  The  effect  of  target  suppression  with 
log  video  is  discussed  in  the  section  on  target  resolution  (Table  4). 

Nonparametric  Detectors 

Nonparametric  detectors  obtain  CFAR  usually  by  ranking  the  test  sample  with  the  reference 
cells  [30,31],  Specifically,  one  orders  the  samples  from  the  smallest  to  the  largest  and  replaces  the 
smallest  with  the  rank  0,  the  next  smallest  with  rank  1 ,  . . .  ,  and  the  largest  with  rank  n  -  1 .  Under 
the  hypothesis  that  all  the  samples  are  independent  samples  from  an  unknown  density  function,  the 
test  sample  has  equal  probability  of  taking  on  any  value.  For  instance,  referring  to  the  ranker  in 
Fig.  20,  the  test  cell  is  compared  to  15  of  its  neighbors.  Since  in  the  set  of  16  samples  the  test  sam¬ 
ple  has  equal  probability  of  being  the  smallest  sample  (or  equivalently  any  other  rank),  the  probability 
that  the  test  sample  takes  on  values  0,  1,  . . .  ,  15  is  1/16.  A  simple  rank  detector  is  constructed  by 
comparing  the  rank  to  a  threshold  K,  and  then  generating  a  1  if  the  rank  is  equal  or  larger,  a  0  if  the 
rank  is  smaller.  The  0s  and  Is  are  summed  in  a  moving  window.  This  detector  incurs  a  CFAR  loss 
of  about  2  dB  but  achieves  a  fixed  Pt„  for  any  unknown  noise  density  as  long  as  the  time  samples  are 
independent.  This  detector  was  incorporated  into  a  postprocessor  used  in  conjunction  with  the 
Federal  Aviation  Administration’s  ASR  radars.  The  major  shortcoming  of  this  detector  is  that  it  is 
fairly  susceptible  to  target  suppression  (e.g.,  if  a  large  target  is  in  the  reference  cells,  the  test  cell 
cannot  receive  the  highest  ranks'. 

If  the  time  samples  arc  correlated,  the  rank  detector  will  not  yield  CFAR.  A  modified  rank 
detector,  called  the  modified  generalized  sign  test  (MGST)  [26],  maintains  a  low  Pf»  and  is  shown  in 
Fig.  21.  This  detector  can  be  divided  into  three  parts:  a  ranker,  an  integrator  (in  this  case  a  two- 
pole  filter),  and  a  threshold  (decision  process).  A  target  is  declared  when  the  integrated  output 
exceeds  two  thresholds.  The  first  threshold  is  fixed  (equals  /t  +  T\/K ,  Fig.  19)  and  yields 
Pfa  -  10  6  when  the  reference  cells  are  independent  and  identically  distributed.  The  second  thres¬ 
hold  is  adaptive  and  maintains  a  low  Pfo  when  the  reference  samples  are  correlated.  The  device  esti¬ 
mates  the  standard  deviation  of  the  correlated  samples  with  the  mean  deviate  estimator,  where 
extraneous  targets  in  the  reference  cells  have  been  excluded  from  the  estimate  by  use  of  a  preliminary 
threshold  T2. 

The  rank  and  MGST  detectors  are  basically  two-sample  detectors.  They  decide  that  a  target  is 
present  if  the  ranks  of  the  test  cell  are  significantly  greater  than  the  ranks  of  the  reference  cells.  Tar¬ 
get  suppression  occurs  at  all  interfaces  (e.g.,  land,  sea),  where  the  homogeneity  assumption  is 
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Fig  20  —  Rank  detector:  output  of  a  comparator  C 
is  either  a  zero  or  one  From  G.V.  Trunk  17) 
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violated.  However,  some  tests  exist,  such  as  Spearman  Rho  and  Kendall  Tau  tests,  that  depend  only 
on  the  test  cell,  these  tests  work  on  the  principal  that  as  the  antenna  beam  sweeps  over  a  point  tar¬ 
get.  the  signal  return  increases  and  then  decreases.  Thus,  for  the  test  cell,  the  ranks  should  follow 
the  pattern  ot  tirst  increasing  and  then  decreasing.  Although  these  detectors  do  not  require  reference 
cells  and  hence  have  the  useful  property  of  not  requiring  homogeneity,  they  are  not  generally  used 
because  of  the  large  CFAR  loss  taken  for  moderate  sample  sizes.  For  instance,  the  CFAR  losses  are 
approximately  10  dB  for  16  pulses  on  target  and  6  dB  for  32  pulses  on  target  [32], 

A  basic  disadvantage  of  all  nonparametric  detectors  is  that  one  loses  amplitude  information, 
which  can  be  a  very  important  discriminant  between  target  and  clutter  [33],  For  example,  a  large 
return  (a  >  lOOOtrr)  in  a  clutter  area  is  probably  just  clutter  breakthrough.  See  “contact  entry 
logic"  in  the  section  entitled  Track-Whilc-Scan  System. 

Clutter  Mapping 

A  clutter  map  uses  adaptive  thresholding  where  the  threshold  is  calculated  from  the  return  in  the 
test  cell  on  previous  scans  rather  than  from  the  surrounding  reference  cells  on  the  same  scan.  This 
technique  has  the  advantage  that  for  essentially  stationary  environments  (e.g.,  land-based  radar  against 
ground  clutter),  the  radar  has  inter  clutter  visibility — it  can  see  between  large  clutter  returns.  At  the 
Lincoln  Laboratory  [34]  they  very  effectively  used  a  clutter  map  for  the  zero  doppler  filter  in  their 
Moving  Target  Detection  (MTD).  The  decision  threshold  T  for  the  ith  pulse  is 

T=ASi_ ,  (23) 

where 


S,  =  K  Sj  _  |  +  X,.,  (24) 

Sj  is  the  background  level,  X,  is  the  ith  pulse,  K  is  the  feedback  value  that  determines  the  map  time 
constant,  and  A  is  the  constant  that  determines  the  false  alarm  rate.  In  the  MTD,  used  for  ASR 
application,  K  is  7/8,  which  effectively  averages  the  last  eight  scans.  Several  questions  concerning 
the  use  of  clutter  maps  remain  unanswered:  (1)  if  frequency  agility  is  used,  will  there  be  a  need  to 
use  a  map  for  each  frequency  or  can  an  average  clutter  map  be  used;  (2)  when  there  is  slow  moving 
clutter,  such  as  rain  and  chaff,  how  should  one  update  the  map  accounting  for  the  clutter  that  moves 
from  one  cell  to  another;  and  (3)  can  clutter  maps  be  used  effectively  on  slowly  moving  platforms 
such  as  ships. 

Target  Resolution 

In  automatic  detection  systems,  a  single,  large  target  will  probably  be  detected  many  times;  e.g., 
in  adjacent  range  cells,  azimuth  beams,  and  elevation  beams.  Therefore,  automatic  detection  systems 
have  algorithms  for  merging  the  individual  detections  into  a  single  centroided  detection.  However, 
most  algorithms  have  been  designed  so  that  they  will  rarely  ever  split  a  single  target  into  two  targets, 
thus  resulting  in  poor  range  resolution  capability.  A  merging  algorithm  commonly  used  [35]  is  the 
adjacent-detection  merging  algorithm,  which  decides  whether  a  new  detection  is  adjacent  to  any  of  the 
previously  determined  set  of  adjacent  detection:-.  If  the  new  detection  is  adjacent  to  any  detection  in 
the  set  of  adjacent  detections,  it  is  added  to  the  set.  Two  detections  are  adjacent  if  two  of  their  three 
parameters  (range,  azimuth,  and  elevation)  are  the  same  and  the  other  parameter  differs  by  the  resolu¬ 
tion  element:  range  cell  A R,  azimuth  beamwidth  0,  or  elevation  bandwidth  7.  A  simulation  was  run 
to  compare  the  resolving  capability  of  some  common  detection  procedures  shown  in  Fig.  22;  and  the 
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Eic.  ;2  —  V'iria(ions  of  basic  detector  T  is  the  threshold,  m  is  either  the  mean  of  all  the  reference  cells  or  is  the  smaller 
(the  minimum!  of  the  means  of  the  reference  cells  at  only  the  right  and  of  the  reference  cells  at  only  the  left.  F  is  a  fixed 
number,  and  n  •>.  the  variance,  (copyright  1978.  IEEE).  From  G  V  Trunk  (35) 


results  are  shown  in  Table  3.  As  demonstrated,  the  log  detectors  have  much  better  performance  than 
the  detectors  using  linear  video.  The  effect  of  target  suppression  on  log  video  can  be  seen  in  Table  4. 
One  notes  an  improved  performance  when  one  calculates  the  threshold  using  only  the  half  of  the 
reference  with  the  minimum  mean  value.  Figure  23  shows  the  resolution  capability  of  the  log  detec¬ 
tor  that  uses  only  the  reference  cells  with  the  minimum  mean.  The  probability  of  resolving  two  equal 
amplitude  targets  does  not  rise  above  0.9  until  they  are  separated  in  range  by  2.5  pulsewidths. 

Assuming  that  the  target  is  small  with  respect  to  the  pulscwidth  and  that  the  pulse  shape  is 
known,  the  resolution  capability  can  be  improved  by  fitting  the  pulse  shape  to  the  data  and  comparing 
the  residue  square  error  to  a  threshold  [36],  If  only  one  target  is  present,  the  residue  should  only  be 
noise  and  hence  should  be  small.  If  two  or  more  targets  are  present,  the  residue  will  contain  signal 
and  should  be  large.  The  results  of  resolving  two  targets  with  S/N  =  20  dB  are  shown  in  Fig.  24. 
These  targets  can  be  resolved  at  a  resolution  probability  of  0.9  at  separations  that  vary  between  1/4 
and  3.4  of  a  pulsewidth  depending  on  the  relative  phase  difference  between  the  two  targets.  Further¬ 
more,  this  result  can  be  improved  further  by  processing  multiple  pulses. 

Detection  Summary 

When  only  2  to  4  samples  (pulses)  are  available,  a  binary  integrator  should  be  used  to  avoid 
detections  from  interference.  When  a  moderate  number  of  pulses  (5  to  16)  are  available,  a  binary' 
integrator,  rank  detector,  or  moving  window  integrator  should  be  used.  If  the  number  pulses  is  large 
(greater  than  20),  a  batch  processor  or  two-pole  filter  should  be  used.  If  the  samples  are  independent, 
a  one  parameter  threshold  can  be  used.  If  the  samples  are  dependent,  one  can  either  use  a  two 
parameter  threshold  or  adapt  a  one-parameter  threshold  on  a  sector  basis.  These  rules  should  only 
serve  as  a  general  guideline.  It  is  highly  recommended  that  the  radar  video  from  the  environment  of 
interest  be  collected  and  analyzed  and  that  various  detection  processes  be  simulated  on  a  computer  and 
tested  against  recorded  data. 
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Table  3  —  Number  of  Targets  Detected  when  Two  Targets  are  Separated  by  1.5,  2.0,  2.5,  3.0, 
or  3.5  Range-Cells  and  the  Third  Target  is  7.0  Cells  from  the  First  Target.  From  Trunk  [35) 


TARGET  SEPARATION  IN  RANGE  CELL  DIMENSIONS  Ah 

Fig.  23  --  Resolution  capability  of  log  detector  that  uses  reference  cells 
with  minimum  mean,  (copyright  1978,  IEEE).  From  G.V.  Trunk  [35] 


TARGET  SEPARATION  IN  PUlSEW'DThS 

Fig.  24  —  Probability  of  resolution  as  a  function  s.f  range  separation:  coherent  ad  hoc 
procedure;  sampling  rate  A R  =  1.5  samples  per  pulsewidth;  target  strengths — 
nonfluctuating,  A ,  =  A2  =  20  dB;  phase  differences  =  0°,  45°,  90°,  135°,  and  180°, 
(copyright  1984,  IEEE).  From  G.V.  Trunk  [36[ 
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AUTOMATIC  TRACKING 

Track-while-scan  systems  are  tracking  systems  for  surveillance  radars  whose  nominal  scan  time 
(revisit  time)  is  from  4  to  12  s.  If  the  probability  of  detection  ( PD )  per  scan  is  high,  if  accurate  tar¬ 
get  location  measurements  are  made,  if  the  target  density  is  low,  and  if  there  are  only  a  few  false 
alarms,  the  design  of  the  correlation  logic  and  tracking  filter  is  straightforward.  However  in  a  realis¬ 
tic  radar  environment  these  assumptions  are  seldom  valid,  and  the  design  of  the  automatic  tracking 
system  is  complicated.  In  actual  situations  one  encounters  target  fades  (changes  in  signal  strength 
caused  by  multipath  propagation,  blind  speeds,  and  atmospheric  conditions),  false  alarms  (caused  by 
noise,  clutter,  interference,  and  jamming),  and  poor  radar  parameter  estimates  (caused  by  noise, 
unstabilized  antennas,  unresolved  targets,  target  splits,  multipath  propagation,  and  propagation 
effects).  The  tracking  system  must  deal  with  all  these  problems. 

A  generic  track-while-scan  will  be  considered  first.  This  will  be  followed  by  a  discussion  of 
maximum  likelihood,  multiple  scan  approaches  to  automatic  tracking. 

Track- While-Scan  System 

A  general  outline  of  a  track-while-scan  system  is  presented.  Next,  since  a  major  portion  of  any 
tracking  system  must  deal  with  manipulating  a  large  amount  of  data  efficiently,  a  typical  file  system 
and  pointer  system  are  described.  Then  the  contact  entry  logic,  coordinate  systems,  tracking  filter, 
maneuver-following  logic,  track  initiation,  and  correlation  logic  are  discussed  in  detail. 

System  Organization 

Almost  all  track-while-scan  systems  operate  on  an  azimuthal  sector  basis  that  provides  basic  sys¬ 
tem  timing.  Figure  25  shows  a  typical  series  of  operations.  For  instance,  if  the  radar  has  reported 
all  the  detections  in  sector  1 1  and  is  now  in  sector  12,  the  tracking  program  would  start  by  correlating 
(try  ing  to  associate)  the  clutter  points  (stationary  tracks)  in  s  ctor  10  with  detections  in  sectors  9.  10, 
1 1.  Those  detections  that  are  associated  with  clutter  points  a.e  deleted  (are  not  used  for  further  corre¬ 
lations)  from  the  detection  file  and  are  used  to  update  the  clutter  points.  Next,  firm  tracks  in  sector  8 
are  correlated  with  detections  in  sectors  7,  8,  and  9.  By  this  lime  all  clutter  points  have  been 
removed  from  sectors  9  and  below  .  Those  detections  that  are  associated  with  firm  tracks  are  deleted 
from  the  detection  file  and  are  used  to  update  the  appropriate  track. 
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I'Mialh,  some  provision  is  made  lor  jj ivinii  prelerence  to  firm  tracks  (instead  of  tentative  tracks) 
m  the  correlation  process.  My  performing  the  eorrclation  process  two  sectors  behind  firm-track  corre¬ 
lations  d  ig.  25).  it  is  impossible  lor  tentative  tracks  to  steal  detections  belonging  to  firm  tracks.  In 
some  other  tracking  systems  the  correlation  lor  firm  and  tentative  tracks  is  performed  in  the  same  sec¬ 
tor;  however,  the  generalized  distance  /)  (see  Fq.  (49))  between  tracks  and  detections  is  incremented 
by  A/)  if  the  track  is  tentative. 

finally,  detections  that  are  not  associated  with  either  clutter  points  or  tracks  are  used  for  initia¬ 
tion  ol  new  tracks.  The  most  common  initiation  procedure  is  to  initiate  a  tentative  track.  Later  the 
tentative  track  is  dropped  or  else  made  a  firm  track  or  clutter  point  depending  on  its  velocity.  An 
alternate  approach  1 37 ]  is  to  establish  both  a  clutter  point  and  a  tentative  track.  If  the  detection  came 
from  a  stationary  target,  the  clutter  point  will  be  updated  and  the  tentative  track  will  eventually  be 
dropped.  However,  if  the  detection  came  from  a  moving  target,  the  tentative  track  will  be  made  firm 
and  the  clutter  point  will  be  dropped.  This  method  requires  less  computer  computation  time  when 
most  of  the  detections  are  clutter  residues. 

File  System 

When  a  track  is  established  in  the  software  of  the  computer  a  track  number  is  assigned  to  it. 
All  parameters  associated  with  a  given  track  are  referred  to  by  this  track  number.  Typical  track 
parameters  are  smoothed  and  predicted  position  and  velocity,  time  of  last  update,  track  quality, 
covariances  matrices  if  a  Kalman  type  filter  that  is  being  used,  and  track  history  (i.e.,  the  last  n 
detections).  Each  track  number  is  also  assigned  to  a  sector  so  that  the  correlation  process  can  be  per¬ 
formed  efficiently  [381.  In  addition  to  the  track  file,  a  clutter  file  is  maintained.  A  clutter  number  is 
assigned  to  each  stationary  or  very  slowly  moving  target.  All  parameters  associated  with  a  clutter 
point  are  referred  to  by  this  clutter  number.  Again,  each  clutter  number  is  assigned  to  a  sector  in 
azimuth  for  efficient  correlation. 

The  basic  files  and  pointer  systems  are  described  below  in  a  FORTRAN  format.  Note  that 
several  higher  level  languages  such  as  PASCAL  have  pointer  systems  that  permit  efficient  implemen¬ 
tation. 

Track  and  Clutter  Number  Files— Only  the  operation  of  the  track  number  file  is  described,  since 
the  operation  of  clutter  number  file  is  identical.  The  development  in  Ref.  38  is  used.  The  track 
number  parameters  are  listed  in  Table  5. 


Table  5  —  Track  Number  Parameters.  From  Cantrell  et  al.  [38] 


Parameters  Description 

NT  Track  number 

DROPT  I  (obtain)  or  0  (drop)  a  track  number  NT 

FULLT  Number  of  available  track  numbers 

NEXTT  Next  track  number  available 

LASTT  Last  track  number  not  being  used 

LISTT(M)  File  whose  M  locations  correspond  to  track  numbers 

M  Maximum  number  of  tracks 
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The  track  number  file  is  initialized  by  setting  LISTT(I)  =  1  +  1  for  1  =  1  through  M. 
LISTT(M)  is  set  equal  to  zero  (denoting  the  last  available  track  number  in  the  file),  NEXTT  =  1  (the 
next  available  track  number),  LASTT  =  M  (the  last  track  number  not  being  used),  and  FULLT  = 
M  -  1  (indication  that  M  -  1  track  numbers  are  available).  Figure  26  shows  a  flowchart  of  the 
operation. 


Fig.  26  -•  Flowchart  for  track  number  file  From 
B.H.  Cantrell,  G.V.  Trunk,  and  J.D.  Wilson  (38) 


When  a  new  track  number  is  requested,  DROPT  is  set  equal  to  one,  and  the  system  checks  to 
sec  if  FULLT  is  zero.  If  FULLT  is  not  equal  to  zero  the  routine  is  called.  Since  DROPT  =  1,  the 
new  track  is  assigned  the  next  available  track  number;  i.e.,  NT  =  NEXTT.  The  next  available  track 
number  in  the  list  is  found,  and  NEXTT  is  set  equal  to  LISTT(NT).  FULLT  is  decremented,  thus 
indicating  that  one  less  track  number  is  available.  Finally,  LISTT(NT)  is  set  equal  to  512  (a  number 
larger  than  the  number  of  possible  tracks).  This  is  not  necessary  but  sometimes  helps  in  debugging 
the  program. 

When  DROPT  ~  0,  a  track  number  NT  is  dropped  by  setting  the  last  available  track  number 
LISTT(LASTT)  equal  to  the  track  number  NT.  LISTT(NT)  is  set  equal  to  zero  to  denote  the  last 
track  number,  and  LASTT  is  then  set  equal  to  the  track  number  being  dropped,  LASTT  •-  NT.  The 
parameter  FULLT  is  incremented,  indicating  that  one  more  track  number  is  available.  Thus,  the 
(rack  (and  clutter)  number  files  maintain  a  linkage  from  one  number  to  the  next,  thereby  eliminating 
searching  techniques. 

Track  Number  Assignment  to  Azimuth  Sector  Files— The  azimuth-range  plane  is  usually  separated 
into  64  or  128  equal  azimuth  sectors.  After  a  track  is  updated  or  initiated,  the  predicted  position  of 
the  target  is  cheeked  to  see  which  sector  it  occupies,  and  the  track  is  assigned  to  this  sector,  If  the 
track  is  dropped  or  moves  to  a  new  sector,  it  is  dropped  out  of  the  sector  in  which  it  was  previously 
located.  The  parameters  associated  with  the  track  sector  files  arc  listed  in  Table  6.  Only  the  assign¬ 
ment  of  track  numbers  to  azimuth  sectors  is  described,  since  the  clutter  number  assignment  is  identi¬ 
cal.  The  TBX(I)  file  contains  the  first  track  number  in  sector  I.  If  TBX(I)  =  0.  no  tracks  arc  in  sec¬ 
tor  I.  The  IDT(M)  file  has  storage  locations  corresponding  to  each  of  the  possible  M  track  numbers. 
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Table  b  Track  Sector  1'ilcs.  From  Cantrell  et  al.  |38| 


Parameter 

I  Jeseription 

TBX(I) 

IDT(M) 

First  track  numbet  in  sector  I  (a  subscript  of 
array  ID'Ti 

Hach  location  corresponds  to  a  track  number, 
and  the  location  contains  the  next  track  number 
in  sector  I  or  a  zero. 

The  first  track  number  in  sector  1  is  obtained  from  FIRST  =  TBX(I).  The  second  track  number  in 
the  sector  is  obtained  by  NFXT1  =  IDT(FIRST).  The  next  track  number  in  the  sector  is  obtained  by 
NHXT2  =  lDT(NFXTl).  The  process  is  continued  until  a  zero  is  encountered,  indicating  that  no 
more  track  numbers  are  in  the  sector. 

When  a  new  track  is  added  or  a  track  moves  from  one  ector  to  another,  a  track  number  must 
be  added  to  the  new  sector.  The  flowchart  tor  achieving  this  is  shown  in  Fig.  27.  The  first  track 
number  in  the  sector  is  stored,  the  track  number  NT  being  added  is  made  the  first  track  number  in  the 
sector,  and  the  track  number  NT  in  the  IDT(.NT)  file  is  made  equal  to  the  original  first  track  in  the 
sector.  This  procedure  is  essentially  a  push-down  stack  that  pushes  the  older  track  numbers  further 
down  in  the  file. 


Fig.  27  —  Adding  a  track 
number  NT  to  sector  1.  From 
B.H.  Cantrell.  C-.V.  Trunk,  and 
J.D.  Wilson  (38| 


When  a  track  is  dropped  or  moves  out  of  the  sector,  the  track  number  must  be  removed  from 
the  sector.  Figure  28  shows  the  flow  diagram  for  this.  First,  it  is  determined  whether  the  first  track 
number  in  the  sector  TBX(I)  is  the  one  being  dropped.  If  it  is,  the  first  track  number  in  the  sector  is 
set  equal  to  the  second  track  number  in  the  sector,  and  the  location  in  the  IDT  file  corresponding  to 
the  track  number  NT  being  dropped  is  set  to  zero.  NT  is  set  equal  to  the  track  number  in  the  file  fol¬ 
lowing  the  one  just  dropped,  so  that  we  now  have  the  next  available  track  number.  If  the  track 
number  being  dropped  is  not  the  first  one  in  the  file,  then  the  push-down  stack  IDT(NL)  is  searched 
sequentially  until  the  track  number  is  found.  The  variable  IDT(NL)  containing  NT  as  the  next  track  is 
replaced  by  the  next  track  number  following  NT,  and  the  variable  in  the  IDT  file  corresponding  to 
NT  is  set  equal  to  zero.  Again,  NT  is  set  equal  to  the  track  number  in  the  file  after  the  one  being 
dropped. 
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Fig.  28  —  Dropping  a  track  number  NR  from  sector  I.  From 
B.H.  Cantrell,  G.V.  Trunk,  and  J.D.  Wilson  [38] 


Contact  Entry1  Logic 

Not  all  the  detections  declared  by  the  automatic  detector  are  used  in  the  tracking  process;  rather, 
many  of  the  detections  (contacts)  arc  filtered  out  in  software  using  a  process  called  “contact  entry 
logic"  [33 1 .  The  basic  idea  is  to  use  the  target  amplitude  in  connection  with  the  type  of  environment 
the  detection  comes  from  to  eliminate  some  detections.  The  first  step  is  for  the  operator  to  identify 
the  environment  in  various  regions.  Types  of  environments  that  can  be  identified  are  land  clutter, 
rain  clutter,  sea  clutter,  and  interference.  Histograms  of  detections  in  various  regions  are  developed 
and  Fig.  29  shows  an  example  for  rain  clutter.  In  general,  the  amplitude  of  rain  clutter  detections  is 
lower  than  that  of  targets.  Consequently  a  threshold  can  be  set  to  inhibit  low-amplitude  detections. 
This  threshold  can  be  controlled  by  the  local  detection  density;  raised  in  high-detection  densities  and 
lowered  in  low-detection  densities.  In  no  circumstances  would  a  detection  be  inhibited  if  it  fell  within 
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SIGNAL/NOISE  (dB) 

Fig.  29  —  Effectiveness  of  signal/noise  ratio  test  in  rain  clutter.  From 
From  W.G.  Bath,  L.A.  Biddson,  S.F.  Haase,  and  E.C.  Wetzlar  [33] 


a  track-gate  (i.e.,  gate  centered  about  the  predicted  position  of  a  firm  track).  In  land  clutter  areas, 
high-amplilude  detections  would  be  inhibited. 

Coordinate  Systems 

The  target  location  measured  by  the  radar  is  in  spherical  coordinates:  range,  azimuth,  elevation, 
and  possibly  range  rate  (doppler).  Thus,  it  may  seem  natural  to  perform  tracking  in  spherical  coordi¬ 
nates.  However,  this  causes  difficulties  since  the  motion  of  the  constant-velocity  targets  (straight 
lines)  will  cause  acceleration  terms  in  all  coordinates.  A  simple  solution  to  this  problem  is  to  track  in 
a  Cartesian  coordinate  system.  Although  it  may  appear  that  tracking  in  Cartesian  coordinates  will 
seriously  degrade  the  accurate  range  track,  it  has  been  shown  [39]  that  the  inherent  accuracy  is  main¬ 
tained.  Another  approach  [40]  noted  that  since  maneuvering  targets  cause  a  large  range  error  but  a 
rather  insignificant  azimuth  error,  a  target-oriented  Cartesian  coordinate  system  could  be  used. 
Specifically,  the  x  axis  is  taken  along  the  azimuth  direction  of  the  target  and  the  v  axis  is  taken  in  the 
cross-range  direction.  Finally,  no  matter  what  coordinate  system  is  used  for  tracking,  the  correlation 
of  detections  with  tracks  must  be  performed  in  spherical  coordinates. 
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Tracking  Filters 

The  simplest  tracking  filter  is  the  a-j3  filter  described  by 

v(  (k  )  =  xp  (k )  4-  ce[jcm(A: )  -  *,,(*)], 


Klk)  =  -  l)  +  flxa(k)  -  xp(k)]/T ,  (26) 

xp(k  +  1)  =  xjk)  +  Vs(k )  T ,  (27) 

where  xjk)  is  the  smoothed  position.  1^, (A: )  is  the  smoothed  velocity.  xp(k)  is  the  predicted  position, 
xm(k)  is  the  measured  position,  T  is  the  scanning  period  (time  between  detections),  and  a  and  /3  are 
the  system  gains. 

The  minimal  mean-square  error  (MSE)  filter  for  performing  the  tracking  when  the  equation  of 
motion  is  known  is  the  Kalman  filter,  first  discussed  by  Kalman  |41]  and  later  by  Kalman  and  Bucy 
J42] .  The  Kalman  filter  is  the  most  commonly  used  filter  in  radar  and  is  a  recursive  filter  that 
minimizes  the  MSE.  The  state  equation  in  xy  coordinates  for  a  constant-velocity  target  is 


X(t  +  1)  =  0(f)  X(t)  +  T(0  A(t) , 


where 


r  i 

x(n 

1  TOO 

x(t) 

0  100 

V(()  ■  - 

oo  i  r 

_  vH)  J 

00  0  1 

T2/ 2 

0 

ax(t) 

T 

0 

0 

T2/ 2 

,  and  A  (t)  = 

0 

T  J 

a,  ( t ) 

no  = 


with  Xi'n  being  the  state  vector  at  time  t  consisting  of  position  and  velocity  components  x(t),  x (/), 
y(0  and  y(i):  t  +  1  being  the  next  observation  time;  T  being  the  time  between  observations;  and 
«,(/)  and  a%U)  being  random  accelerations  with  covariance  matrix  Q(t).  The  observation  equation  is 
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(32) 

with  Yu  i  being  the  measurement  at  time  t  consisting  of  positions  .v,„(r)  and  y,„ (/ )  and  V(t)  being  a 
/cro-mcan  noise  whose  covariance  matrix  is  R(t). 
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The  problem  is  solved  recursively  by  first  assuming  the  problem  is  solved  at  time  t  -  1. 
Specifically  ,  it  is  assumed  that  the  best  estimate  X(t  -  1  i  t  -  1)  at  time  r  -  1  and  its  error  covari¬ 
ance  matrix  P(t  —  1  |  r  —  1)  are  known,  where  the  caret  in  the  expression  of  the  form  X(t  |s) 
signifies  an  estimate  and  the  overall  expression  signifies  X(t)  is  being  estimated  with  observations  up 
to  y(^).  The  six  steps  involved  in  the  recursive  algorithm  are: 

•  calculate  the  one-step  prediction 

X(r|r  -  1)  =  $>(/  -  l)X(r  -  l|r  -  1);  (33) 

•  calculate  the  covariance  matrix  for  the  one-step  prediction 

P(t  |r  -  1)  =  <*>(/  -  \)  P(t  -  1 1 r  -  l)4>r(/  -  1)  +  T(r  -  1  )Q(t  -  l)rr(r  -  1);  (34) 

•  calculate  the  ,  redicted  observation 

Y(t\t  -  1)  =  M{t)X(t\t  -  1);  (35) 


•  calculate  the  filter  gain 

A(r)  =  P(r  jr  -  1)  MT(t)  [ M(t)P(t\t  -  \)MT{t)  +  R( OP1;  (36) 


•  calculate  the  new  smoothed  estimate 

X(t  |  r)  =  X(t  \t  -  1)  +  A(r)  [K(r)  -  Y(t  \  t  -  1));  and  (37) 

•  calculate  the  new  covariance  matrix 

P(t  |0  =  (/  -  A(OM(0]  Pit  \t  ~  1).  (38) 


In  summary  ,  starting  with  an  estimate  X{t  -  l)rr  —  1)  and  its  covariance  matrix  P(t  —  1 1  r  —  1 ) 
after  a  new  observation  Y(t)  is  received  and  the  six  quantities  in  the  recursive  algorithm  are  calcu¬ 
lated,  a  new  estimate  X(r  1 1)  and  its  covariance  matrix  P(t  \  t)  are  obtained. 


It  can  be  shown  [43]  that  for  a  zero  random  acceleration,  Q{t)  s  0  and  a  constant  measurement 
covariance  matrix  R(t)  =  R ,  the  a-Q  filter  can  be  made  equivalent  to  the  Kalman  filter  by  setting 


and 


2(2 k  -  1) 
k(k  +  1) 


(39) 


on  the  A.  th  scan. 


d  = 


6 

k  (k  +  1 ) 


(40) 


Thus  with  time,  tc  and  ,3  approach  zero,  and  heavy  smoothing  is  applied  to  the  new  samples. 
Usually  it  is  worthwhile  to  bound  a  and  d  from  zero  by  assuming  a  random  acceleration  Q{t )  s*  0 
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corresponding  to  approximately  a  1-g  maneuver.  The  previously  stated  Kalman  filter  method  is 
optimal  (with  respect  to  MSE)  for  straight-line  tracks  but  must  be  modified  to  enable  the  filter  to  fol¬ 
low  target  maneuvers. 

Maneuver- Following  Logic 

Benedict  and  Bordner  [44]  noted  that  in  track-while-scan  systems,  there  is  a  conflicting  require¬ 
ment  between  good  tracking  noise  reduction  (implying  small  «  and  &)  and  good  maneuver-following 
capability  (implying  large  a  and  (3).  Although  some  compromise  is  always  required,  the  smoothing 
-quations  should  he  constructed  to  give  the  best  compromise  for  a  desired  tracking  noise  reduction. 
Benedict  and  Bordner  defined  a  measure  of  transient-following  capability  and  showed  that  a  and  i3 
should  be  related  by 


Thus,  an  (a,  d)  pair  satisfying  Eq.  (41)  can  be  chosen  so  that  the  tracking  filter  will  follow  a  speci¬ 
fied  g-turn.  Cantrell  [45]  calculated  the  probability  that  a  target  detection  will  fall  within  a  correla¬ 
tion  region  centered  at  the  predicted  target  position  when  the  target  is  doing  a  specified  g-turn.  He 
suggested  using  the  (a.tf)  pair  satisfying  Eq.  (41)  that  yields  the  smallest  correlation  region.  How¬ 
ever,  if  high  g-turns  must  be  followed,  the  noise  performance  is  poor. 

An  alternate  approach  uses  a  turn  detector  (Fig.  30)  that  consists  of  two  cot  relation  regions.  The 
inner  nonmaneuvering  gale  is  usually  set  so  that  the  probability  of  the  target  detection  being  in  the 
gate  is  greater  than  0.99  when  the  target  is  doing  a  1-g  maneuver.  Then,  if  the  detection  is  in  the 
nonmaneuvering  correlation  region,  the  fiber  operates  as  usual,  the  filter  gains  being  reduced  accord¬ 
ing  to  Eq.  (36)  or  Eqs.  (39)  and  (40).  When  the  target  detection  falls  outside  the  inner  gate  but 
within  the  maneuver  gate,  a  maneuver  is  declared  and  the  filter  bandwidth  is  increased  (a  and  0  are 
increased)— Quigley  and  Holmes  [46]  increase  the  bandwidth  by  lowering  the  value  of  k  in  Eqs.  (39) 
and  (40).  To  avoid  the  problem  of  a  target  fading  causing  a  missed  detection  and  a  false  alarm 
appearing  in  the  large  maneuver  gate,  the  track  should  be  bifurcated  when  a  maneuver  is  declared. 
That  is,  two  tracks  are  generated:  the  continuation  of  the  old  track  but  not  updated  with  the  new 
detection  and  a  new  maneuvering  track  with  the  new  detection  and  increased  filter  bandwidth.  The 
next  scan  is  used  to  resolve  the  ambiguity  and  remove  one  of  the  tracks.  Although  the  turn  detector 
is  the  most  common  approach  for  maneuver-following,  other  solutions  are  to  adjust  the  bandwidth  as 
a  function  of  the  measurement  error  [37]  or  to  use  the  Kalman  filter  with  a  realistic  target- 
maneuvering  model  [47], 


SMOOTHED 

POSITION 


MANEUVERING  GATE 


NONMANEUVERING 

GATE 


Pig.  TO  —  Maneuver  and  nonnianeuvcr  gales  eenlercd  at  the 
target's  predicted  position  Prom  G.V  Trunk  |7| 
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Specifically,  Cantrell  et  al.  [37]  suggested  that  the  a-0  filter  (described  by  Eqs.  (25)  through 
(27))  be  made  adaptive  by  adjusting  a  and  0  by 


1  -X+J 

a  =  1  —  e 

(42) 

and 

0  = 

1  +  c~nu'J  -  2e~^J  cos  (oj0Wl  -  r2), 

(43) 

in  which 

u0  =  0.5  \p  \(k)/p2{k)  |  , 

(44) 

where 

Pi(*)  = 

e  ~u°Tp \(k  -  1)  +  (1  -  e~UbT)(  (k)e(k  -  1), 

(45) 

p2(k) 

=  p2(k  -  1)  +  (1  -  e~“J)  e(*)e(*). 

(46) 

f  is  the  damping  coefficient 

(nominally  0.7),  T  is  the  time  since  the  last  update, 

and  Wj,  are 

weighting  constants,  and  t(k)  is  the  error  between  the  measured  and  predicted  positions  on  the  Arth 
update  The  basic  principle  of  the  filter  is  that  p{(k)  is  an  estimate  of  the  covariance  of  successive 
errors  and  p2(k)  is  an  estimate  of  the  error  variance.  When  the  target  trajectory  is  a  straight  line, 
p{(k)  approaches  0,  since  the  expected  value  of  e(k)  is  0.  Thus  u0  approaches  0,  and  the  filter  per¬ 
forms  heavy  smoothing.  When  the  target  turns,  p^k)  grows,  since  the  error  e(A )  will  have  a  bias. 
Thus  «0  grows,  and  the  filter  can  follow  the  target  maneuver. 

Singer  [47]  suggested  using  the  Kalman  filter  with  a  realistic  target-maneuvering  model.  He 
assumed  that  the  target  was  moving  at  a  constant  velocity  but  was  perturbed  by  a  random  accelera¬ 
tion.  The  target  acceleration  is  correlated  in  time;  and  it  was  assumed  that  the  covariance  of  the 
correlation  was 


r(r)  =  E\a(l)a{t  +  r)|  =  a^e  ^  |  r  ]  , 


(47) 


where  a(r)  is  the  target  acceleration  at  time  /,  o*  is  the  variance  of  the  target  acceleration,  and  7  is 
the  reciprocal  of  the  maneuver  time  constant.  The  density  function  for  target  acceleration  consists  of 
delta  functions  at  ±<4max,  a  delta  function  at  0  with  probability  P0,  and  a  uniform  density  between 
and  /4nm.  For  this  density 
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For  this  target  motion.  Singer  then  calculated  the  state  transition  matrix  <t>(t)  and  the  covariance 
matrix  Q(t ).  thereby  specifying  the  Kalman-filter  solution.  He  generated  curves  that  give  the  steady- 
state  performance  of  the  filter  for  various  data  rates,  single-look  measurement  accuracies,  encounter 
geometry,  and  class  ot  maneuvering  targets. 


G. V  TRUNK 


Track  Initiation 

Detections  that  do  not  correlate  with  clutter  points  or  tracks  are  used  to  initiate  new  tracks.  If 
the  detection  does  not  contain  doppler  information,  the  new  detection  is  usually  used  as  the  predicted 
position  (in  some  military  systems,  one  assumes  a  radially  inbound  velocity);  and  a  large  correlation 
region  must  be  used  for  the  next  observation.  The  correlation  region  must  be  large  enough  to  capture 
the  next  detection  of  the  target  that  has  the  maximum  velocity  of  interest.  Since  the  probability  of 
obtaining  false  alarms  in  the  large  coliection  region  is  sometimes  large,  one  should  generally  disre¬ 
gard  the  initial  detection  if  no  correlation  is  obtained  on  the  second  scan.  Also  one  should  not  declare 
a  track  firm  until  at  least  a  third  detection  (falling  within  a  smaller  correlation  region)  is  obtained.  A 
common  track  initiation  criteria  is  four  out  of  five,  although  one  may  require  only  three  detections  out 
of  five  in  regions  with  low  false  alarm  and  target  density  rates.  The  possible  exceptions  for  using 
only  two  detections  are  when  doppler  information  is  available  (so  that  a  small  correlation  region  can 
be  used  immediately  and  the  range  rate  can  be  used  as  an  additional  correlation  parameter)  or  for 
“pop-up"  targets  (i.e.,  targets  that  suddenly  appear  at  a  close  range)  in  a  military  situation. 

An  alternative  track  initiation  logic  [40]  uses  a  sequential  hypothesis-testing  scheme.  When  a 
correlation  is  made  on  the  ith  scan,  A,  is  added  to  a  score  function;  when  a  correlation  opportunity  is 
missed,  yt  is  subtracted  from  the  score  function.  The  increments  are  a  function  of  the  state  of  the 
tracking  system,  the  closeness  of  the  association,  the  number  of  false  alarms  in  the  region,  the  a  priori 
probability  of  targets,  and  the  probability  of  detection.  When  the  score  function  exceeds  a  particular 
value,  the  track  becomes  firm.  Although  this  method  will  inhibit  false  tracks  in  dense  detection 
environments,  it  will  not  necessarily  establish  the  correct  tracks. 

To  initiate  tracks  in  a  dense  detection  environment,  the  technique  known  as  “retrospective  pro¬ 
cessing"  [48]  uses  the  detections  over  the  last  several  scans  to  initiate  straight-line  tracks  by  using  a 
collection  of  filters  matched  to  different  velocities.  Figure  31  illustrates  an  example  of  the  process¬ 
ing;  here  one  looks  for  surface  targets  in  the  presence  of  sea  spikes.  It  is  clear  from  this  figure  that 
detections  1,  4,  6,  10,  12,  and  14  form  a  track.  Even  though  the  false  alarm  rate  per  scan  is  rela¬ 
tively  high  (approximately  10“-\  the  processor  can  be  implemented  with  a  microprocessor. 

Track  Drop 

Firm  tracks  that  are  not  updated  in  several  scans  corresponding  to  40  to  60  s  usually  are 
dropped.  In  some  systems,  before  a  track  is  dropped  the  track  symbol  is  blinked  thus  indicating  a 
track  is  about  to  be  dropped.  This  gives  an  operator  the  chance  to  update  the  track  with  a  detection 
and  keep  it  from  being  dropped. 

Correlation  Logic 

Several  procedures  are  used  to  associate  detections  with  tracks.  Of  special  interest  are  the  con¬ 
flicting  situations  of  multiple  tracks  competing  for  a  single  detection  or  of  multiple  detections  lying 
within  a  track’s  correlation  gate. 

To  limit  the  number  of  detections  that  can  update  a  track,  correlation  gates  are  used.  A  detec¬ 
tion  never  will  update  a  track  unless  it  lies  within  the  correlation  gate  that  is  centered  at  the  track’s 
predicted  position.  The  correlation  gate  should  be  defined  in  r  —  6  coordinates,  regardless  of  what 
coordinate  system  is  used  for  tracking.  Furthermore,  the  gate’s  size  should  be  a  function  of  the  mea¬ 
surement  accuracy  and  prediction  error  (specified  by  Eq.  (34))  so  that  the  probability  of  the  correct 
detection  lying  within  the  gate  is  high  (at  least  0.99).  In  some  tracking  systems  (49),  the  location  of 
the  correlation  gate  is  fed  back  to  the  automatic  detector;  and  the  detection  threshold  is  lowered  in  the 
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gate  to  increase  the  probability  of  detection  PD.  The  gate  also  disables  the  “contact  entry’’  logic 
described  in  this  section. 

When  several  detections  are  within  the  correlation  regions,  the  usual  and  simplest  solution  is  to 
associate  the  closest  detection  with  the  track.  Specifically,  the  measure  of  closeness  is  the  statistical 
distance 


Dz  = 


(r»  “  rmY 


-  + 


o; 


(0p  ~  Am)2 


(49) 


where  ( rp ,  6p)  is  the  predicted  position,  (rm,  0m)  is  the  measured  position,  a ;  is  the  variance  of 
ri>  ~  rm-  an^  's  variance  of  6p  -  d„, .  These  variances  are  a  by-product  of  the  Kalman  filter. 
Since  the  prediction  variance  is  proportional  to  the  measurement  variance,  and  0;  are  sometimes 
replaced  by  the  measurement  variances.  Statistical  distance  rather  than  Euclidean  distance  must  be 
used  because  the  range  accuracy  is  usually  much  better  than  the  azimuth  accuracy. 
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Figure  32  illustrates  problems  associated  with  multiple  detections  and  tracks;  two  detections  are 
within  gate  1,  three  detections  arc  within  gate  2.  and  one  detection  is  within  gate  3.  Table  7  lists  all 
detections  within  the  tracking  gate,  and  the  detections  are  entered  in  the  order  of  their  statistical  dis¬ 
tance  from  the  track.  Tentatively,  the  closest  detection  is  associated  with  each  track,  and  then  the 
tentative  associations  are  examined  to  remove  detections  that  are  used  more  than  once.  Detection  8. 
which  is  associated  with  tracks  1  and  2,  is  paired  with  the  closest  track  (track  1  in  this  case);  then  all 
other  tracks  are  reexamined  to  eliminate  all  associations  with  detection  8.  Detection  7  is  associated 
with  tracks  2  and  3;  the  conflict  is  resolved  by  pairing  detection  7  with  track  2.  When  the  other  asso¬ 
ciation  with  detection  7  is  eliminated,  track  3  has  no  associations  with  it  and  consequently  will  not  be 
updated  on  this  scan.  Thus  track  1  is  updated  by  detection  8,  track  2  is  updated  by  detection  7,  and 
track  3  is  not  updated. 
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An  alternate  strategy  is  to  always  pair  a  detection  with  a  track  if  there  is  only  one  correlation 
with  a  track.  As  before,  ambiguities  are  removed  by  using  the  smallest  statistical  distance.  Thus, 
track  3  in  the  example  is  updated  by  detection  7,  track  1  is  updated  by  detection  8.  and  track  2  is 
updated  by  detection  9.  This  procedure  yields  better  results  when  PD  is  close  to  one  and  the  proba¬ 
bility  of  false  alarm  is  very  low  . 
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Table  7  —  Association  Table  for  the  Hxamplc 
Shown  in  Fig.  32.  From  Trunk  |7|. 
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The  optimal  strategy  that  maximizes  the  probability  of  correct  correlation  is  a  joint  maximum 
likelihood  approach.  This  involves  examining  all  possible  combinations  and  selecting  the  combination 
that  is  most  probable  in  a  statistical  sense.  This  procedure  requires  knowing  the  probability  of  detec¬ 
tion  and  probability  of  false  alarms.  A  discussion  of  the  joint  maximum  likelihood  approach  can  be 
found  in  the  section  entitled  Maximum  Likelihood  Approaches. 

Singer  and  Sea  [50]  were  two  of  the  first  people  to  recognize  and  characterize  the  interaction 
between  the  correlation  and  track  update  functions.  Specifically,  three  distinct  situations  can  occur: 
the  track  is  not  updated,  the  track  is  updated  with  the  correct  return,  and  the  track  is  updated  with  an 
incorrect  return.  They  generalized  the  tracking  filter’s  error  covariance  equations  to  account  for  the  a 
priori  probability  of  incorrect  returns  being  correlated  with  the  track.  This  permits  the  analytical 
evaluation  of  tracking  accuracy  in  a  multitarget  environment  that  produces  false  correlations.  Further¬ 
more,  by  using  the  generalized  tracking  error  covariance  equation,  they  optimized  the  lilter  gain 
matrix  that  yielded  a  new  minimum-error  tracking  filter  for  multitarget  environments.  Also,  they 
generated  a  suboptimal  fixed-memory  version  of  this  filter  to  reduce  computation  and  memory 
requirements. 

A  subsequent  paper  by  Singer  et  al.  [51)  uses  a  posteriori  correlation  statistics  based  on  all 
reports  in  the  vicinity  of  the  track.  The  mathematical,  structure  is  similar  to  that  of  the  Kalman  filter. 
The  estimation  error  is  denoted  by  X(t  [/')  =  X(t)  —  X(t  \t ')  and  has  mean  and  covariance 
matrices  denoted  by  b(t  |r')  and  Pit  |  / ').  It  is  assumed  that  nk  detections  fall  within  the  correlation 
gate  on  scan  k .  Included  in  the  number  nk  are  extraneous  reports  whose  number  obeys  a  Poisson  dis¬ 
tribution  and  whose  positions  are  uniformly  distributed  within  the  gate.  The  smooth  estimate  is  given 
by 


XU  |0  =  X(t  j/  -  !)  +  <4(0,  (50) 

where  A  (t)  is  chosen  to  minimize  the  noncentral  second  moment  of  the  filter  estimation  error. 

The  problem  is  solved  by  using  track  histories.  A  track  history  a  at  scan  k  is  defined  by  select¬ 
ing,  for  each  scan,  a  detection  or  a  miss  (a  miss  corresponds  to  the  hypothesis  that  none  of  the  reports 
belong  to  the  track).  The  number  of  such  track  histories  is 

Uk)  =  £  (I  +  «,).  (51) 

i  - 1 

A "undated  with  each  history  n  is  the  probability  />„(/)  that  the  history  a  is  the  correct  one,  given 
observations  through  time  i  (scan  k).  The  terms  (t\t  -  1)  ai  j  P'lt{t  |f  -  D  arc  <hc  bias  and 
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covariance  of  the  estimation  error  X(t  ]  t  -  1),  given  observations  through  time  ?  —  1  and  given  that 
track  history  a'  at  time  t  -  I  is  the  (only)  correct  one.  Recursive  equations  are  obtained  for  pa,  ba, 
and  then  as  shown,  the  optimal  correction  vector  is  given  by 

A(t)  =  £  pa(t)ba(t\t  -  1).  (52) 

a  — 

This  solution  not  only  minimizes  the  mean-square  error  but  also  is  an  unbiased  estimate. 

The  optimal  a  posteriori  filter  requires  a  growing  memory;  therefore  several  suboptimal  filters 
were  investigated  by  Singer  et  al.  [51].  The  first  suboptimal  filter  considers  only  the  las*  .V  scans; 
track  histories  that  are  identical  for  the  last  N  scans  are  merged.  The  second  suboptimal  filter  only 
considers  the  L  nearest  neighbors  in  the  correlation  gate;  essentially  the  gate  size  is  changed  to  limit 
the  number  of  reports  to  L.  The  last  method  uses  both  techniques:  it  considers  only  the  last  ;V  scans 
and  restricts  the  number  of  reports  on  any  scan  to  L  . 

Figures  33  and  34  show  the  results  of  simulations  that  compare  the  optimal  and  suboptimal  a 
posteriori  filters,  the  optimal  and  suboptimal  a  priori  filters,  and  the  Kalman  filter.  In  Fig.  33,  the 
filter  variance  normalized  by  the  theoretical  (perfect-correlation)  Kalman-filter  variance  is  plotted  for 
several  filters.  As  a  class  the  a  posteriori  filters  provide  better  performance  than  the  other  filters. 
However,  for  high  density  of  false  reports  (4(3 a\  -  0.1),  the  a  posteriori  filter  is  30  times  worse  than 
predicted  by  the  standard  Kalman-filter  approach.  Thus  the  standard  approach  should  never  be  used 
in  dense-target  (or  false-target'  environments.  Figure  34  gives  the  probability  of  making  a  false 
correlation.  Again  the  a  posteriori  filters  provide  the  best  performance. 

If  one  considers  only  the  last  scan  (.V  -  0  case),  one  obtains  the  probabilistic  data  association 
filter  (PDAF)  [52,53].  In  essence,  each  detection  is  used  to  update  a  Kalman  filter  and  the  final  esti¬ 
mate  is  the  weighted  sum  of  all  the  estimates  of  each  detection  where  the  weight  is  the  probability  that 
cadi  detection  is  the  proper  update  This  procedure  does  not  require  a  growing  memory;  it  requires 
slightly  higher  computational  requirements  than  the  ordinary  Kalman  filter  and  solves  the  correlation 
problem  by  using  all  detections  with  the  track  gate  to  update  the  filter.  As  originally  formulated, 
track  i n it  ation  and  termination  were  not  considered.  However,  Colcgrove  et  al.  [54,55]  have  pro¬ 
posed  soli  lions  to  these  problems.  Furthermore,  they  have  implemented  the  PDAF  on  an  OTH  radar. 

Maximum  Likelihood  Approaches 

The  previously  discussed  procedures  updated  a  single  fack  with  cither  zero,  one,  or  multiple 
detections;  however,  each  track  was  operated  on  individually;  i.e.,  each  track  was  considered  by 
itself.  Better  performance  can  be  obtained  by  using  all  the  detections,  on  all  the  previous  scans  by 
employing  a  maximum  likelihood  approach.  This  better  performance  is  obtained  at  the  cost  of  an 
enormous  computational  requirement. 

Sittler  |56]  formulated  the  maximum  likelihood  approach  to  tracking.  In  his  approach  he 
assumed  that:  (a)  objects  appeared  in  the  surveillance  region  according  to  a  Poisson  model  and  were 
uniformly  distributed  in  >pacc.  (b)  the  persistence  of  objects  is  independent  with  a  duration  that  fol¬ 
lows  an  exponential  density,  (e)  the  track  positions  have  a  known  density,  (d)  the  detection  model  is  a 
Poisson  process,  (c)  the  measurement  errors  have  a  known  density  (later  assumed  to  be  Gaussian), 
and  (f)  the  false  alarms  follow  a  Poisson  process. 

Stein  and  Blackman's  |57]  approach  is  similar  to  that  of  Sittler  [56].  It  differs  in  the  following 
ways:  (a)  u  sample  data  system  is  assumed  and  the  detection  process  is  specified  by  the  detection 
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33  -  Simulated  and  calculated  variances  in  filtered  position  errors  for  optimal 
priori  tracking  filters,  (copyright  1974,  IHKE).  From  R.A.  Singer,  R.G.  Sea, 
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probability,  (b)  the  development  eliminates  the  concept  of  undetected  targets  and  of  distinguishing 
between  true  targets  generating  a  single  detection  and  false  alarms,  (c)  it  uses  realistic  target 
maneuver  models  [45],  and  (d)  it  uses  the  Kalman  filter  [41,42],  Later,  Morefield  [58]  developed  a 
similar  approach  suitable  for  implementation.  It  is  based  upon  likelihood  functions  and  converts  the 
association  of  detections  to  tracks  into  an  integer  programming  problem.  Finally,  Trunk  and  Wilson 
[59]  considered  the  effects  of  resolution  in  the  track  formation  process. 

Since  the  general  approach  of  all  these  maximum  likelihood  procedures  is  similar,  the  more 
comprehensive  approach  [59]  will  be  discussed  in  detail.  The  maximum  likelihood  method  calculates 
the  likelihood  that  a  set  of  tracks  correctly  represents  a  given  set  of  detections.  The  statement  that  a 
set  of  tracks  represents  a  given  set  of  detections  means  that  each  detection  is  either  declared  a  false 
alarm  or  assignee  to  one  or  more  tracks  and  only  one  detection  is  assigned  to  a  track  on  any  scan.  In 
the  calculation  of  the  likelihood,  the  probability  of  detection,  the  probability  of  false  alarm,  the  mea¬ 
surement  error  characteristics,  and  the  probability  of  target  resolution  are  all  taken  into  account.  For 
simplicity,  the  likelihood  is  only  written  in  terms  of  the  range  measurement.  Reference  59  indicates 
how  the  angle  measurements  are  included. 

The  probability  of  obtaining  Nj  detections  of  the  j'th  target  in  Ns  scans  is  binominally  distri¬ 
buted.  Then  the  probability  of  obtaining  the  detections  associated  with  a  set  of  1 V,  tracks  is  just  the 
product  of  the  binomial  probabilities  associated  with  each  track,  i.e., 


PDN,(  1  -  Pd)n,~n‘- 


(53) 


The  probability  of  obtaining  a  specified  number  of  false  alarms  in  a  range  interval  R,  on  each  of 
the  Ns  scans  is  given  by  the  multinomial  probability 


n,\  n  ^ 

l  - 1 


where  ,\D  is  the  total  number  of  detections,  M{  is  the  number  of  scans  that  contain  exactly  i  false 
alarms,  and  P,  is  the  Poisson  probability  of  obtaining  /  false  alarms  on  any  scan,  i.e.. 

P,  --  (\R,V  exp  (-XT?,)//!  (55) 


where  X  is  the  false  alarm  rate  per  unit  length.  Finally,  the  probability  of  a  false  alarm  occurring  at  a 
particular  range  is  given  by  the  uniform  density.  Consequently  the  likelihood  of  a  specified  set  of 
false  alarms  is  obtained  by  multiplying  Eq.  (54)  by 

.v, 

(\/R,)  £  i  M,.  (56) 

/  =  i 


The  likelihood  of  obtaining  a  given  set  ot  measurements  on  the  yth  track  is 


V.-  i 

n  f\XjU).  Xlruc(/ )]/(/•’,  )2  Gwcrf  ~2,u 


(57) 
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where  a2  is  the  measurement  error  variance,  Xj(i)  is  the  range  of  the  associated  detection  on  the  tth 
scan,  Xtrue(/')  is  the  true  range  of  the  track  on  the  ith  scan,  and  the  function /(•,  )  is  defined  by 


1 ,  if  there  is  no  detection  associated 
with  track  on  /  th  scan  or  if  the 
associated  detection  is  also  associated 
with  any  other  track 

exP(-(.r  -  y)2/2 o2],  otherwise. 


(58) 


The  rationale  behind  Eq.  (57)  is  as  follows.  The  detections  occurring  on  the  earliest  and  latest  scan, 
which  may  not  be  the  first  and  the  last  scans,  can  occur  anywhere  in  the  range  interval  this  in 
essence  establishes  an  a  priori  velocity  distribution.  Hence,  the  factor  (\/R,)2  reflects  the  uniform 
distribution  of  the  first  and  last  detections.  The  other  detections  are  Gaussian  distributed  with  mean 
value  equal  to  the  true  position  and  a  variance  of  a2.  Since  the  true  positions  are  unknown,  the  obvi¬ 
ous  approach  is  to  replace  the  true  positions  by  the  minimum  mean  square  estimate  |X;(£)].  To  avoid 
a  bias  caused  by  the  estimation  procedure,  it  was  shown  [59]  that  the  initial  and  final  error  measure¬ 
ments  should  be  included  in  Eq.  (57).  Consequently,  one  should  use 

n  f[XJ(i),XJ{i)/{Rl)2(hTo2f  -2>/2]  (59) 

1  =  I 


instead  of  Eq.  (57). 

When  several  targets  are  close  to  one  another,  merging  algorithms  may  yield  a  single  detection 
from  the  closely  spaced  targets  [35]  The  probability  of  obtaining  a  single,  unresolved,  detection  Xk 
from  \'k  closely  spaced  targets  is  calculated  by  first  ordering  the  predicted  positions  so  that 

4  s  4  s  •••  *  4 


where  for  notational  convenience  the  scan  identifier  has  been  dropped.  Letting  Dm  =  X ^  , 

the  probability  of  not  resolving  the  Nk  targets  is  given  by 


p r (Xfc )  =  n  P(Dm) , 

m  =  2 


(60) 


where  the  probability  of  not  resolving  two  targets  separated  by  distance  D  is  given  in  Ref.  35  by 


P(D)  = 


1, 

(2.6 R  -D)/0.9R , 

0, 


D  <  MR 
\.1R  <  D  <  2.6 R 
D  <  2.6 R 


(61) 
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where  R  is  the  3  dB  pulsewidth  (range  cell  dimension).  It  was  shown  that  the  density  of  the  position 
of  Xk  can  be  approximated  by 

Pp(Xk)  ~  [1/(2t(t2)1/2]  PE(Xk),  (62) 

where 

PE(Xk)  =  [exp  (  — €i/2<r2]/max[l,(XIw  -  Xl)/(2ira2)u2} ,  (63) 

i 

and 

e*  =  max  (0,  Xk  -  X,v ,  X,t  -  Xk)  (64) 

i 


is  the  distance  from  Xk  to  the  nearest  predicted  position  if  Xk  lies  outside  of  the  inteival  defined  by 
the  predicted  positions;  otherwise  (k  is  zero.  Since  the  factor  (2 iro2)172  in  Eq.  (62)  is  already 
included  in  Eq.  (59),  the  probability  associated  with  unresolved  detections  is  found  by  multiplying  Eq. 
(60)  and  Eq.  (63). 

The  final  term  in  the  likelihood  includes  the  a  priori  probability  of  having  Np  tracks  in  the  local 
region.  Since  in  some  situations  this  can  be  the  result  of  a  deterministic  decision  (e.g.,  a  wing  of  air¬ 
craft  in  a  military  situation),  the  a  priori  probabilities  will  be  assumed  to  be  equal  rather  than  the 
commonly  used  Poisson  arrival  model  [56,57].  In  applying  the  method,  any  a  priori  probability 
appropriate  to  the  problem  of  interest  can  be  used. 

In  terms  of  the  previous  expressions,  the  likelihood  of  an  N p  track  combination  is  given  by 


s, 

L(Nt)  = 

n 

.  j= 1 

J  j 

PdN'(  1  -  Pd)N,~N' 

•  jo/*/) 


'*0  " 

E  i  M, 

ns  !  n 
/■  =  1 

0=i 

‘  1 

,vr  V, 


(XR,)iMl  e  ,!(/!)"' 

11 


n  n/[(x/iu;(i)]/(/?/)2(2fJ2) 


{N  —  2)/2 


;  =  i  h=l 


iJ 


n  PR(Xk)pE(Xk) 

*  =  i 


(65) 


where  ,YU  is  the  number  of  unresolved  detections. 

Figure  35  shows  an  example  of  this  procedure  applied  to  recorded  data  of  a  flight  of  five  air¬ 
craft.  The  five  tracks  identified  are  ( 1 ,  M ,  1 ,  1 ,  M ,  1 ),  (M ,  1,2,2,  1 ,  2),  (2,  2,  2,  M ,  1 ,  2),  (3,  3, 
3,  3.  2.  3),  and  (4,  4,  4,  4,  3,  4);  each  track  is  represented  by  a  sextuple  with  the  number  specifying 
a  detection  number  and  an  M  specifying  a  miss.  We  note  that  tracks  2  and  ^  have  common  detec¬ 
tions  on  scans  3.  5.  and  6.  Using  Eq.  (65),  this  sequence  of  tracks  was  over  100  times  more  likely 
than  the  best  four  track  combination. 
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hood  tracks 


MILTISENSOR  INTEGRATION 

Mr\isensor  integration  will  be  divided  into  three  topics:  (1)  colocated  radar  integration,  (2) 
multisite  radar  integration,  and  (3)  integration  of  data  from  unlike  sensors  (e  g.,  radar  and  IFF).  In 
this  section,  radars  will  be  considered  to  be  colocated  if  it  is  not  necessary  to  take  into  account  that 
the  radars  are  located  a:  different  sites. 


There  are  several  methods  for  integrating  data  from  multiple  radars  into  a  single  system  track 
file.  The  type  of  radar  integration  that  should  be  used  is  a  function  of  the  radar’s  performance,  the 
environment,  and  whether  or  not  the  radars  are  colocated.  The  several  integration  methods  that  have 
been  used  in  various  systems  are: 

•  Track  Selection 

Generate  a  track  with  each  radar  and  choose  one  of  the  tracks  as  the  system  track. 

•  Average  Track 

Generate  a  track  with  each  radar  and  weight  according  to  the  Kalman  filter’s  covariance 
matrices  the  individual  tracks  to  form  a  system  track. 

•  Augmented  Track 

Generate  a  track  with  each  radar,  choose  one  of  the  tracks  as  the  system  track,  and  use 
selected  detections  from  the  other  radars  to  update  the  system  track. 

•  Deteclion-to-Track 

Use  all  radar  detections  to  update  the  system  track;  tracks  may  or  may  not  be  initiated  using 
all  detections  from  all  radars. 

Theoretically,  the  detection-to-track  method  of  integration  yields  the  best  tracks  because  all  the 
available  information  is  used.  However,  the  dete'tions  must  be  weighted  properly  and  care  must  be 
taken  so  that  bad  data  do  not  corrupt  good  data. 
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There  are  many  advantages  in  radar  integration.  Probably  the  most  important  is  that  it  provides 
a  common  surveillance  picture  to  all  users  so  that  decisions  can  be  made  more  effectively.  Radar 
integration  will  also  improve  track  continuity  and  tracking  on  maneuvering  targets  because  of  the 
higher  effective  data  rate.  Improvement  in  track  initiation  times  is  a  function  of  the  target  trajectory. 
For  instance,  long-range  targets  are  usually  detected  by  only  one  radar  so  that  little  or  no  improve¬ 
ment  in  initiation  time  is  achieved.  However,  there  could  be  an  appreciable  reduction  in  the  initiation 
time  for  pop-up  targets.  Finally,  the  general  tracking  performance  is  improved  in  an  ECM  environ¬ 
ment  because  of  the  integration  of  radars  in  different  frequency  bands  located  at  different  positions— 
providing  both  spatial  and  frequency  diversity. 

The  main  advantage  of  integrating  different  sensors  with  radar  is  that  it  provides  classification 
and  or  identification  information  on  radar  tracks.  In  general,  the  other  sensors  do  not  provide  posi¬ 
tion  data  of  accuracy  comparable  to  the  radar  data.  In  addition  to  providing  classification  informa¬ 
tion,  the  sensors  can  alert  each  other  of  conditions  that  can  cause  the  mode  of  operation  to  be 
changed.  For  instance,  a  strong  direction  finding  (DF)  bearing  strobe  that  cannot  be  correlated  with 
any  radar  track,  may  cause  the  radar  to  use  burn  through,  to  lower  its  detection  thresholds,  or  to 
change  its  initiation  criterion  in  the  sector  containing  the  DF  bearing  strobe. 

The  next  three  sections  describe  colocated  radar  integration,  multisite  radar  integration,  and 
unlike  sensor  integration. 

Colocated  Radar  Integration 

The  fundamental  work  [37,60,61]  on  colocated  radar  integration  has  been  performed  for  the 
U.S.  Navy;  typically,  with  a  naval  ship  that  has  two  or  three  radars  within  several  hundred  feet  of 
one  another.  Although  various  radar  integration  techniques  have  been  investigated,  the  one  imple¬ 
mented  on  most  naval  ships  is  the  detection-to-track  integration  philosophy.  Starting  with  the  DDG- 
2  15  class  of  ships  [60,61]  the  U.S.  Navy  started  placing  Integrated  Automatic  Detection  and  Track¬ 
ing  tIADT)  systems  aboard  various  classes  of  U.S.  Navy  ships. 

Figure  36  shows  the  typical  functions  of  an  IADT  employing  the  detection-to-track  integration 
philosophy.  An  automatic  detector  that  performs  a  thresholding  operation  to  control  the  false  alarm 
rate  is  associated  with  each  radar.  After  detections  (contacts)  have  been  declared,  some  further  tests 
(contact  entry  logic)  are  applied  to  see  whether  the  detections  should  be  made  available  to  the  tracking 
system  or  discarded.  For  instance,  a  detection  falling  within  a  tracking  gate  is  never  discarded,  while 
a  detection  within  an  area  that  has  been  declared  to  be  land  clutter  is  discarded  if  its  amplitude  is  too 
large.  (For  a  fuller  description  of  the  contact  entry  logic  see  the  section  entitled  Contact  Entry 
Logic.)  As  detections  become  available  to  the  tracking  system,  in  turn  they  are  associated  with  the 
stationary-track  filter  (clutter  map),  the  track  file,  and  the  saved  detection  file  (previously  uncorrelated 
detections).  Thus,  a  detection  from  one  radar  can  be  correlated  with  a  saved  detection  from  another 
radar,  resulting  in  a  track  entry  sooner  than  possible  from  any  individual  radar. 

The  tracking  algorithms  for  a  multiple  radar  tracking  system  are  quite  similar  to  those  for  a  sin¬ 
gle  radar  tracking  system.  Consequently,  we  will  discuss  only  the  difference. 

First,  a  time  is  associated  with  each  detection;  then,  the  detection  with  the  oldest  detection  time 
is  subjected  to  the  correlation  process.  In  a  single  radar  tracking  system,  time  is  not  required  since 
all  the  timing  can  take  place  in  terms  of  the  radar  s.  in  time. 

The  second  difference  is  that  the  tracking  fiber  updates  the  target  tracks  with  detections  whose 
accuracy  varies  with  the  radar  and  that  arrive  “randomly”  in  time.  The  optimal  tracking  filter  (with 
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Fiji  76  —  Detection-to-track  processing,  (copyright  1977,  IEEE). 
From,  P  G  Casner,  Jr  ,  and  R.J  Prengaman  |60) 


respect  to  the  minimum  mean  square  error)  is  the  Kalman  filter  [41,42],  Consequently,  either  a  Kal¬ 
man  filter  or  an  approximation  of  it  should  be  used. 

Another  problem  unique  to  the  multiradar  tracking  system  is  the  azimuthal  and  range  alignment 
between  radars  This  problem  is  solved  by  adaptively  aligning  all  the  radars  to  the  3-D  radar  that  is 
arbitrarily  selected  to  be  the  accurate  radar.  This  is  accomplished  by  examining  the  radar  tracks  that 
are  being  updated  by  multiple  radars  and  that  meet  other  specified  criteria  (e  g.,  frequency  of  update). 
Then  the  average  error  for  the  tth  radar,  over  all  targets  between  the  measured  and  predicted  posi¬ 
tion,  is  weighted  and  subtracted  from  future  measurements  made  by  the  tth  radar,  over  all  targets. 
This  feedback  arrangement  will  drive  the  bias  errors  toward  zero.  The  residue  bias  error  is  approxi¬ 
mately  110  the  accuracy  of  the  measurement  and  prediction  errors  ]60,61]. 

Multi-Site  Radar  Integration 

The  methods  used  for  exchanging  tracking  information  between  noncolocated  sites  are  a  function 
of  whether  the  sites  are  fixed  or  mobile  and  what  communication  links  are  available. 

The  most  common  multiple  site  tracking  system  philosophy  is  that  one  site  controls  a  radar  track 
(has  reporting  responsibility),  transmits  the  tracking  parameters  over  the  communication  link,  and 
receiving  units  display  the  remote  track.  The  track  receives  only  updates  from  one  site  and  thus  does 
not  benefit  from  the  mutual  support  available  from  other  radars  ,n  different  frequencies  and  spatial 
location.  This  philosophy  is  usually  used  in  any  multisite  system  when  the  communication  link  has 
limited  capacity. 

Multisitc  integration  with  large  gridlock  and  sensor  misalignment  errors  is  a  difficult  problem. 
II  the  targets  from  the  multiple  sites  can  be  correlated,  then  one  can  use  the  Kalman  filtering  tech¬ 
niques  to  estimate  the  gridlock  and  misalignment  errors  (62).  However,  if  the  gridlock  and  misalign¬ 
ment  errors  are  large,  one  cannot  perform  the  necessary  correlation.  Bath  [63]  solved  this  problem 
by  solving  the  correlation,  gridlock,  and  misalignment  problems  simultaneously.  The  maximum 
errors  m  latitude,  longitude,  and  azimuth  misalignment  must  first  be  defined.  Next,  one  must  define 
error  bins  of  si/c  At  in  latitude.  Ay  in  longitude,  by  10  in  azimuth  and  divide  the  error  space  into 
error  bins  of  this  size.  The  error  between  one  track  from  site  1  and  one  track  from  site  2  is  calcu¬ 
lated.  and  a  "one”  is  entered  into  the  bin  corresponding  to  the  calculated  error.  If  there  are  n  tracks 
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at  site  1  and  m  tracks  at  site  2,  m  x  n  errors  are  calculated  and  m  x  n  ones  are  entered  into  the 
error  bins.  For  tracks  that  are  different,  the  ones  will  be  randomly  distributed  throughout  the  error 
space.  However,  for  tracks  that  are  the  same,  all  the  ones  will  appear  in  the  same  error  cell.  An 
example  of  this  two-dimensional  cross-correlation  function  is  shown  in  Fig.  37  where  the  azimuth 
error  is  zero.  The  spike  not  only  gives  the  gridlock  biases  but  also  identifies  the  proper  correlation  of 
tracks  from  multiple  sites.  Kalman  filter  techniques  could  be  used  with  the  correlated  track  pairs  to 
obtain  more  accurate  estimates  of  the  biases. 


Fig  37  —  Two-dimensional  cross-correlation  function  (for  two  track  bases 
with  x  and  v  biases  only).  From  W'.G.  Bath  [63] 


The  netted  radar  program  NRP  (64]  is  a  netted  radar  system  for  ground-target  surveillance 
developed  at  the  Lincoln  Laboratory  under  the  sponsorship  of  the  U.S.  Army  and  the  Defense 
Advanced  Research  Projects  Agency  (DARPA).  The  basic  idea  is  to  associate  automatic  detection 
and  tracking  systems  (similar  to  the  MTD  system)  with  each  radar  and  then  transmit  tracks  to  a  cen¬ 
tral  location  (Radar  Control  Center)  where  the  tracks  are  integrated  into  a  single  track  file.  Tracks 
aie  transmitted  rather  than  detections  because  the  bandwidth  required  to  transfer  tracks  is  significantly 
lower  than  that  required  to  transfer  detections.  The  bandwidth  used  in  the  NRP  system  is  approxi¬ 
mately  2400  baud.  The  NRP  system  was  tested  at  Fort  Sill,  Oklahoma.  The  automatic  integration  at 
the  Radar  Control  Center  provided  continuity  in  track  updating  and  identity  despite  missing  detections 
caused  by  terrain  masking  and  occasional  misses  caused  by  a  variety  of  sources.  Note  that  gridlock 
was  cbtained  by  surveying  the  radar  sites.  This  is  the  usual  and  preferred  solution  to  gridlock  when 
the  sites  are  fixed  and  there  arc  time  and  resources  to  perform  the  survey. 
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The  communication  link  is  one  of  the  keys  to  a  successful  multisite  integration  system.  This  is 
not  a  problem  for  traffic  control  and  land-based  air  defense  systems  that  can  use  telephone  lines;  how¬ 
ever,  it  can  be  a  problem  for  other  systems.  Usually  communication  links  are  not  adequate  to 
transmit  all  the  useful  data  that  are  available  at  the  sensor.  A  very  effective  alternate  way  of  transmit¬ 
ting  data  is  by  means  of  the  radar.  By  using  the  radar  transmitter  and  antenna,  essentially  jam-proof 
communications  can  be  established  (65). 

Unlike  Sensor  Integration 

A  number  of  sensors  can  be  integrated:  radar,  IFF,  IR,  optical,  and  acoustic.  However, 
although  there  has  been  much  talk  and  several  programs  in  sensor  integration,  there  has  been  very  lit¬ 
tle  progress  made.  Electromagnetic  sensors  such  as  radar,  IFF,  and  strobe  extractors  appear  to  be 
most  easily  integrated. 

IFF  Integration 

The  problem  of  integrating  radar  and  IFF  data  is  less  difficult  than  that  of  integrating  two 
radars.  The  question  of  whether  detections  or  tracks  should  be  integrated  is  a  function  of  the  applica¬ 
tion.  In  a  military  situation,  by  integrating  detections,  one  could  interrogate  the  target  only  a  few 
times,  identify  it,  and  then  associate  it  with  a  radar  track.  From  then  on,  there  would  be  little  need 
for  reinterrogating  the  target.  However,  if  tracks  were  integrated,  targets  would  be  interrogated  many 
times  resulting  in  possibly  revealing  their  position.  However,  in  an  air  traffic  situation,  targets  would 
be  interrogated  every  scan  and  consequently  either  detections  or  tracks  can  be  integrated. 

Radar-DF  Bearing  Strobe  Integration 

Coleman  [66]  and  later  Trunk  and  Wilson  [67,68]  considered  the  problem  of  correlating  radar 
tracks  with  DF  bearing  strobes.  Trunk  and  Wilson  considered  the  problem  of  given  K  DF  angle 
tracks,  each  specified  by  a  different  number  of  DF  measurements,  associate  each  DF  track  with  either 
no  radar  track  or  one  of  m  radar  tracks,  again  each  radar  track  being  specified  by  a  different  number 
of  radar  measurements.  Since  each  target  can  carry  multiple  emitters,  i.e.,  multiple  DF  tracks  can  be 
associated  with  each  radar  track,  each  DF  track  association  can  be  considered  by  itself,  thus  resulting 
in  K  disjoint  association  problems.  Consequently,  an  equivalent  problem  is;  given  a  DF  track  speci¬ 
fied  by  n  DF  bearing  measurements,  associate  the  DF  track  with  no  radar  track  or  one  of  m  radar 
tracks,  the  jlh  radar  track  being  specified  by  m/  radar  measurements.  Thus,  the  problem  reduces  to 
the  multiple  hypothesis  testing  problem. 

DF  measurements  associate  with  no  radar  track 

H,  DF  measurements  associate  with  ;th  radar  track. 

J  J 

Using  a  combination  of  Bayes  and  Neyman  Pearson  procedures  and  assuming  that  the  DF  measure¬ 
ments  errors  are  usually  independent  and  Gaussian  distributed  with  zero  mean  and  constant  variance 
<r  but  with  occasional  outliers  (i.e..  large  errors,  not  described  by  the  Gaussian  density),  they  argued 
that  the  decision  should  be  based  on  the  probability 

P,  =  Probability  (Z  >  d;) ,  (66) 

where  Z  has  a  chi-square  density  with  n)  degrees  of  freedom  and  ds  is  given  by 

n 

dj  =  £  Min[4,  \6eU,)  -  $,-(/,  )]:/<C  1  j  =  D - m  .  (67) 

i  -- 1 
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where  rtj  is  the  number  of  DF  measurements  overlapping  the  time  interval  for  which  the  jih  radar 
track  exists,  0e(t,)  is  the  DF  measurement  at  time  0, (/, )  is  the  predicted  azimuth  of  radar  track  j 
for  time  t, ,  and  the  factor  4  limits  the  square  error  to  4 a2  to  account  for  DF  outliers.  Using  the  larg¬ 
est  P/ s  designated  P ^  and  Pnext  and  thresholds  TL,  TH,  TM,  and  R,  the  following  decision  rule  was 
generated 

Firm  Correlation:  P max  >  TH  and  Pmax  >  Pnext  +  R, 

Tentative  Correlation:  TH  >  /,max  >  TM  and  Pmax  S  Pntx,  +  R ; 

Tentative  Correlation  with  Some  Track :  P  *  tm  but  Pmax  <  Pnext  +  R  , 

Tentatively  Uncorreiated:  TM  >  Pmax  >  TL ;  and 
Firmly  Uncorrelated:  TL  >  P^, 
where  the  decision  means  the  following: 

Firm  Correlation ,  DF  signal  goes  with  the  radar  track  that  has  the  largest  Py  (i.e.,  Pnax) 
Tentative  Correlation,  DF  signal  probably  goes  with  the  radar  track  that  has  the  largest  P;  (i.e.. 


Tentative  Correlation  with  Some  Track,  DF  signal  probably  goes  with  some  radar  track,  but  it 
cannot  determine  which 

Tentatively  Uncorreiated ,  DF  signal  probably  does  not  go  with  any  radar  track 

Firmly  Uncorreiated,  DF  signal  does  not  go  with  any  radar  track. 

The  lower  threshold  TL  determines  the  probability  that  the  correct  radar  track  (i.e.,  the  one 
associated  with  the  DF  signal)  will  be  incorrectly  rejected  from  further  consideration.  This  threshold 
is  set  by  noting  that  the  probability  P;  for  the  correctly  associated  radar  track  is  uniformly  distributed 
between  zero  and  one.  Thus,  if  one  desires  to  keep  the  rejection  rate  for  the  correct  track  below  PR, 
one  can  obtain  this  by  setting  TL  =  PR .  The  threshold  TH  is  set  equal  to  PFa  defined  as  the  proba¬ 
bility  of  falsely  associating  a  radar  track  with  a  DF  signal  when  the  DF  signal  docs  not  belong  with 
the  radar  track.  The  threshold  TH  :  a  r'-.K,;">n  of  the  azimuth  difference  n  between  the  true  (DF) 
position  and  the  radar  track  under  consideration.  The  threshold  TH  was  found  for  the  /x  =  1.0  a  and 
/z  =  1.5  a  by  simulation  techniques,  and  the  results  for  PF  «  0.01  are  shown  in  Fig.  38.  Between 
the  high  and  low  thresholds  there  is  a  tentative  region.  The  middle  threshold  divides  the  “tentative" 
region  into  a  tentative-correlated  region  and  a  tentatively  uncorreiated  region.  The  rationale  in  setting 
the  threshold  is  to  set  the  two  associated  error  probabilities  equal  for  a  particular  separation;  i.e., 
Pr  | P,  <  TSf  '■  correct  match|  =  Pr(Py  >  !  incorrect  match).  The  threshold  TM  was  found  by 

using  simulation  techniques  and  is  shown  in  Fig.  39. 

The  purpose  of  the  probability  margin  R  is  to  ensure  the  selection  of  the  proper  DF-radar  asso¬ 
ciation  (avoiding  false  target  classifications)  when  there  are  two  or  more  radar  tracks  close  to  one 
another.  The  correct  selection  is  reached  by  postponing  a  decision  until  the  two  discriminant  proba¬ 
bilities  differ  by  R.  The  value  for  R  is  found  by  specifying  a  probability  of  an  association  error  Pe 
according  to  P,  -  P«|P;njx  —  Pm</  +  R\<  where  Pniiix  corresponds  to  an  incorrect  association  and 
P,u  u  corresponds  to  the  correct  association.  The  probability  margin  R  is  a  function  of  Pe  and  the 


49 


HIGH  THRESHOLD. 


G.  V.  TRUNK 


NUMBER  OF  SAMPLES 


Fig  38  —  High  threshold  vs  number  of  samples  for  two  different  separations.  The  ,r's  are 
the  simulation  results,  (copyright  1987,  IEEE).  From  G.V.  Trunk  and  J.D.  Wilson  [67] 
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lie  79  -  Middle  threshold  vs  number  of  samples  for  two  different  target 
separations  The  x  s  are  the  simulation  results,  (copyright  1987,  IEEE).  From 
G.V.  Trunk  and  J.D.  Wilson  |67| 
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separation  /<  of  the  radar  tracks.  The  probability  margin  R  was  found  for  n  =  0.25  a,  0.50  a,  and 
1.00  (i  by  using  simulation  techniques,  and  Fig.  40  shows  the  results  for  Pr  -  0.01.  Since  the 
curves  cross  one  another,  we  can  ensure  that  P,.  s  0.01  for  any  /t  by  setting  R  equal  to  the  max¬ 
imum  value  of  any  curve  for  each  value  of  n. 


Fig.  40  —  Probability  margin  vs  number  of  DF  measure¬ 
ments  for  three  different  target  separations.  The  O’s, 
X ’s  and  A’s  are  the  simulation  results  for  =  0.25, 
/i  =  0.5,  and  n  -  1.0,  respectively,  (copyright  1987, 
IEEE).  From  G.V.  Trunk  and  J.D.  Wilson  [67] 


The  algorithm  was  evaluated  by  using  simulations  and  recorded  data.  When  the  radar  tracks  are 
separated  by  several  standard  deviations  of  the  measurement  error,  correct  decisions  are  rapidly 
made.  However,  if  the  radar  tracks  are  close  to  one  another,  errors  are  avoided  by  postponing  the 
decision  until  sufficient  data  are  accumulated.  An  interesting  example  is  shown  in  Figs.  41  through 
43.  Figure  41  shows  the  radar  (azimuth)  measurements  of  the  control  aircraft,  the  radar  measure¬ 
ments  of  four  aircraft  of  opportunity  in  the  vicinity  of  the  control  aircraft,  and  the  DF  measurements 
from  the  radar  on  the  control  aircraft.  The  probability  discriminants,  one  without  limiting,  the  other 
with  limiting,  are  shown  in  Figs.  42  and  43,  respectively.  Initially,  an  aircraft  of  opportunity  has  the 
highest  discriminant  probability;  however,  a  firm  decision  is  not  made  since  Pmax  does  not  exceed 
P next  by  the  probability  margin.  After  the  14th  DF  measurement,  the  emitter  is  firmly  correlated  with 
the  control  aircraft.  However,  at  the  18th  DF  measurement,  a  very  bad  measurement  is  made 
'outlier)  and  the  firm  correlation  is  downgraded  to  a  tentative  correlation.  Figure  43  shows  that,  if 
limiting  is  employed,  the  correct  decision  remains  firm. 

In  a  complex  environment  where  there  arc  many  radar  tracks  and  DF  signal  sources,  it  is  quite 
possible  that  many  DF  signals  will  be  assigned  the  category  that  the  DF  signal  probably  goes  with 
some  radar  track.  To  remove  many  of  these  ambiguities,  multiple-platform  DF  operation  can  be  con¬ 
sidered.  The  extension  of  the  previous  procedures  to  multipie-platform  operation  is  straightforward. 
Specifically,  if  9,, ,(/,)  and  0e2Uk )  are  the  DF  angle  measurements  with  respect  to  platforms  1  and  2, 
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Fig  41  —  Radar  detections  o  and  DF  measurements  collected  on  the 
control  aircraft  The  O's.  ,3's  + 's,  and  x 's  are  radar  detections  on 
four  aircraft  ot  opportunity  in  the  vicinity  of  the  control  aircraft,  (copy¬ 
right  198'.  IEEE)  From  G.V.  Trunk  and  .ED.  Wilson  [67] 
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Fig  43  -•  Probability  discriminant  for  experimental  data:  solid  line  is  discriminant,  and  dashed  lines  are 
thresholds  TM  and  TH.  Limiting  used,  (copyright  1987,  IEEE).  From  G.V.  Trunk,  and  J.D.  Wilson  [67] 

and  0y|(f,)  and  dj2{tk)  are  the  estimated  angular  positions  of  radar  track  j  with  respect  to  platforms  1 
and  2,  the  multiplatform  square  error  is  simply 


dj  =  £  Min  (4,  0r ,(/,•)  -  0j\(tj)\2/o}\ 


+  £  Min  I4-  WeiVk)  -  (68) 

*  =  i 

Then,  the  previously  described  proceduie  can  be  used  with  dj  being  defined  by  Eq.  (68)  instead  of 
Eq.  (67). 
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